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ABSTRACT 

We explore the properties of the H/3 emission hne profile in a large, homogeneous and 
bright sample of N'-^ 470 low redshift quasars extracted from Sloan Digital Sky Survey (DR5). 
We approach the investigation from two complementary directions: composite/median spec- 
tra and a set of line diagnostic measures (asymmetry index, centroid shift and kurtosis) in 
individual quasars. The project is developed and presented in the framework of the so called 
4D Eigenvector 1 Parameter Space, with a focus on its optical dimensions, FWHM(H/3) and 
the relative strength of optical Fell (RFe//=W(FeII4434-4684A)AV(H^)). We reenforce the 
conclusion that not all quasars are alike and spectroscopically they do not distribute randomly 
about an average typical optical spectrum. Our results give further support to the concept of 
two populations A and B (narrower and broader than 4000 km s^^ FWHM(H/5), respectively) 
that emerged in the context of 4DE1 space. The broad H/? profiles in composite spectra of Pop- 
ulation A sources are best described by a Lorentzian and in Population B by a double Gaussian 
model. Moreover, high and low accretion sources (an alternative view of the Population A/B 
concept) not only show significant differences in terms of Black Hole (BH) mass and Edding- 
ton ratio Lhoj/L Edd, but they also show distinct properties in terms of line asymmetry, shift 
and shapes. We finally suggest that a potential refinement of the 4DE1 space can be provided 
by separating two populations of quasars at Ri?e//^0.50 rather than at FWHM(H/3)=4000 
km s~^. Concomitantly, the asymmetry and centroid shift profile measures at 1/4 fractional 
intensity can be reasonable surrogates for the FWHM(H/3) dimension of the current 4DE1. 

Key words: galaxies: active, (galaxies:) quasars: emission hnes, (galaxies:) quasars: general 



1 INTRODUCTION 

1. 1 Current View on the Structure of BLR of AGN 

1.1.1 Observational considerations 

A defining characteristic of Type 1 Active Galactic Nuclei (AGN) 
is the presence of broad emission fines in their optical and UV 
spectra. We cannot resolve tfie AGN central regions at other 
than radio (Very Long Baseline Interferomet ry - milliarcsecond 
resolution) or, more recently, infrared (e.g., iRaban et al.l l2009h 
wavelengths. This leaves the structure, geometry, kinematics 
and nebular physics of the Broad Line Region (BLR) beyond 
our direct reach. It is unlikely that we will spatially resolve the 
central regions of more tfian a few low z AGN in the foreseeable 
future. We are left with two practical methods to probe the 
BLR structure and kinematics: a) reverberation studies (e.g.. 



Koratkar & Gaskelll Il99ll: IPetersonI Il993l : iPeterson & Wandell 
19991: IWandel, P eterson & Ma lkanl 1 19991: (KaspietalJ 1 20001: 



Peterson et al.l 12004: Kaspi et al.l l2007l: iBentz et alj 120091) and 



b) single epoch spectroscopy for large numbers of sources. The 
former requires dedicated spectroscopic monitoring to achieve 
temporal resolution of a relatively small number of sources while 
the latter provides kinematic resolution for multitudes of sources. 
The time lag between fine and continuum flux variations (in 
the former approach) provides an estimate of the BLR size. A 
few dozen (mostly low z) sources have been monitored with 
typical BLR radii in the range from a few light-days to several 



typical isLK radii in trie range trom a tew lignt-aays to several 
light-weeks (e.g., IPetersonI 1 19931 : iKaspi et al.ri2000l : iBentz et al.l 



2009). Spectroscopy of large samples provides a much-needed 
empirical clarification since, without a clear understanding of the 
phenomenology, one can build models filled with misconceptions 
about BLR structure/geometry, kinematics and nebular physics. 



* E-mail: zamfiOOl @crimson. ua.edu 



We have been exploring the phenomenology of quasars 
in the context of a 4D Eigenvector 1 parameter space (4DE1; 
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Sulentic, Marziani & Dultzin-Hacvanl2000al; 


Sulentic et al.l2000bl: 


Marziani et alj 200 iL 2003jbl Sulentic et al. 


2007). which serves 



as a spectroscopic unifier/discriminator for all type 1 AGN. We 
take it as self-evident that such a parameter space is needed: 1) 
to emphasize source differences, 2) to help remove the degener- 
acy between the spectroscopic signatures of orientation and physics 
and 3) to provide a context within which to interpret results of 
multiwavelength studies of individual sources. We adopted four 
principal 4DE1 parameters because they were relatively easy to 
measure in large numbers of sources and because they showed 
large intrinsic dispersion and sometimes major differences between 
subtypes. The principal parameters involve measures of: (i) full 
width at half-maximum of broad H/3 (FWHM(H/3)), (ii) equivalent 
width ratio of optical Fell (A4570A blend) and broad H/3, Kpeii = 
W(FeIIA4570A)AV(H/3), (iii) the soft X-ray photon index {V.oft) 
and (iv) CIVA1549A broad-line profile velocity displacement at 
half-m aximum, c(l/2). Various lines of evidence ( Sulentic et al. 
|2007|) suggested a change of source properties near FWHM(H/3) 
= 4000 km s~^, which motivated the introduction of the Popula- 
tion A-B distinction in order to emphasize the differences. This is 
similar to the radio-quiet (R Q) vs. ra dio loud (RL) distinction, but 
appears to be more effective jSulentic et al. ,2002; see also b elow). 

Previous studies JSulentic et al.l l2003l : 

IZamfir. Sulentic & Marzianil 1200 Sl) indicate that RL sources 
show a much more restricted 4DE1 domain space occupation 
(FWHM(H/3) versus Kf^ii) than the RQ majority (~78% of all 
RL quasars belong to Population B domain and ~62% of RQ 
belong to Population A domain). This was in part the motivation 
for our Population A-B distinction. A majority (« 55-60%) of 
quas ars show broad Balmer lines with FW]-lM(H / 3) < 4000 km 
s~^ ("Sulentic, Marziani & Dultzin-HacyarJ l2000al : ISulentic et~al] 
l2000b. : .Zamfir, Sulentic & Marziani 20080 and are labeled Popula- 
tion A. Population A shows: 1) a scarcit y of RL sources (3-4% at 
most: IZamfir. Sulentic & Marzianill2008l) . 2) strong/moderate Fell 
emission, 3) a soft X-r ay excess (perhaps th e consequence of a 
high accretion rate; e.g., IWang & Netz"eil l20031 4) high-ionization 
broad lines (HIL; e.g., CIV A 1549 A) showing blueshift/asymmetry 
and 5) low ionization broad line profiles (LIL) best described by 
Lorentz fits ('Sulentic et al .l2002l ; lBachev et al.l2004ISulentic et all 
[200 7: Marziani et al. 200^. 

Population B sources: 1) include the large majority of 
the RL sources: ~91% of all FRII and ~62% of all core- 
dominated RL quasars or, alternatively, ~17% of all Population 
B s ources are RL quasars (and this is a conservative estimate; 
see IZamfir. Su lentic & Marziani 2008), 2) show weak/moderate 
Fell emission, 3) do not show a soft X-ray excess, 4) do not show 
HIL blueshift/asymmetry and 5) show LIL Balmer lines best fit 
with d ouble Gaussian mod els ( Sulentic et alj 20021 ; iBachev et al.l 



l2004l : ISulentic et ai] |2007|; iMarziani et al.l |2009|) . Fundamental 
differences between the BLR in Population A and B sources (and 
betwee n RL and RQ q u asars) a ppear to be lower electron density 
(e.g., ICollin-Souffrini 1 19871 ; ICollin-Souffrin. Hameurv & Jolvl 
Il988ah in Population B along with a less flattened 



cloud 



in 

distribution 



(e.g., 'Zheng, Sulentic & Binette| 



ISulentic. Marziani & Dultzin-Hac van 2000a; Collin et al] 12006 
Some of the Population B sources with very broad LIL pro- 



^ Assuming that the true parent population of RL quasars are FRII sources 
there should exist for every core-dominated RL source of a certain radio 
luminosity an even weaker radio power FRII, if the unification scenario is 
valid. 



files show double-peaked structure (e.g., lEracleous & HalpemI 
Il994h . but they are rare. Nonetheless, there is an underlying 
idea that a double-peaked signature is hidden in most profiles. 
iMurrav & Ch iang ( 1997) has emphasized that a disk-l-wind model 
favors the regions of low projected velocity in the disk, which 
leads to a predominance of single-peaked profiles. The existence 
of double-peaked profiles is not a nece ssary condition to p rove 
the reality of a disk BLR geometry (e.g.. |Popovic et ai]|2004l . and 
references therein). 

Valuable insights are provided by studies of emission lines of 
different ionization levels: HIL, e.g., CIV A1549A and LIL, e.g., 
Balmer lines. Fell, Mgll doublet A2800A. The HIL and LIL are 
probably produced in distinct regions (e.g., ICollin -Souffrin et all 
il988h) at least in Population A sources as inferred from the dif- 
ference in the time delayed response to continuum variatio ns (see 
Section 4.1 of lSulentic, Marziani & Dultzin-Hacvanl200d3) . While 
CIV is always broader than H/3 (produced at smaller radii?) in 
Population A sources, the two lines tend to scale m ore closely or, 
in so me cases, H/9 tends to show broader profiles iSulentic et al.l 
l2007h in Population B. 

There has been a great deal of disagreement over the meaning 
of often poor spectroscopic data for often small samples. Making 
progress involves a lot of work obtaining spectra with good S/N and 
resolution in order to assemble representative sample. Ideally one 
would like be able to observe the same line(s) over a large redshift 
range, which raises obvious difficulties that IR spectrosc opy has 
begun to alleviate (e.g., S ulentic et al - 2004, 2006; Marzi ani et al.l 
1 20091 and references therein). Or, even better, observe concomi- 
tantly both HIL and LIL for the same source with awareness that 
the B LR in high and low luminosity sources may be different as 
well (Netzer 2003: Ka spi et alj2007l ; lNetzer & Trakhtenbrojl 20071 ; 
IMarziani e t al. 2009). 

We suggest that the best hope of advancing our understanding 
of the AGN phenomenon lies with better empiricism, which is the 
path to model improvement. The first step to progress requires a 
reasonably large (i.e., representative) quasar sample for which uni- 
form and good quality spectra exist. We argue that progress can- 
not come from studies of indiscriminately averaged quasar spectra, 
whether their S/N is high or low. This statement is false only if all 
quasars are virtually identical, which we know to not be true. 

There is a long series of questions in need of answers. Do 
RQ and RL quasars show the same geometry/kinematics? Do opti- 
cal/UV spectra offer predictive power on the likelihood of a source 
being or becoming RL? How many broad line emitting regions exist 
in a source and is it the same number in all sources? Can we dis- 
tinguish high and low accretors spectroscopically? How are broad 
and narrow emission lines related geometrically and kinemati- 
cally? How much variance in nebular physics exists from source- 
to-source? Do quasars change spectroscopically from low to high 
redshift? 



1. 1.2 Geometrical and Kinematical Inferences 

Current best guesses about geometry and kinematics involves 
the idea that all or much of the broad-line emission arises 
from a flattened cloud distribution (possibly an accretion disk) 
and that kinematics are dominated by Keplerian motions, (e.g., 
|Peterson & Wandel 1999, 2000). There are empirical arguments 
in favor of a flattened geometry, which pertain mostly to radio- 
loud (RL) quasars. An observed anticorrelation between the width 
of the Balmer lines (FWHM) and th e core:lobe radio flux den- 
sity ratio R ( IWills&Brown3ll986l) or equivalent alternatives 
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dWills & Brothertoiilll99i ICorbinlll997h can be explained by ob- 
servation of a flattened BLR at different viewing angles. Similar 
conclusions could be drawn from an analysis of the FWHM as a 
function of radio spectral index a djarvis & McLure 2 006) or the 
Balmer decrement as a function of R ( Jackson & Browneill99i[) . In 
superluminal quasars the viewing angle 9 can be determined with 
good accuracy and there are reported correlations between 6 or 
other beaming indicators and Balmer lines FWHM jRokaki et alj 



Volonteri. S ikora & Lasotal |2007|: Idel R Laeos, Padilla & Cord 
2009 ; Tchekhovs kov. Naravan & McKinneyl |2009|) , the 
structure of the accretion disk (e.g., iBallantvnj |2007| : 



l2003l:ISulenticetal 



5003h that favor an axisymmetric BLR dom- 



inated by rotation. Zamfir. Sul entic & MarzianJ ( 1200 Sl) find that 
FRII (Fanaroff-Rilev njFanar off & Rilev.l974i) quasars show both 
broader H/3 profiles (AFWHM « 2000 km s"^) and weaker Fell 
(AR Fe/i ~ 0.2) than th e core/core-jet RL quasars. 

ICoUin et al. I( l200^ offer arguments that inclination effects are 
apparent only in sources with narrow profiles - what we call Popu- 
lation A. Studies in the 4DE1 context suggest that kinematics and 
geometry are very different for Population B sources. Even if RL 
quasars of Population B show correlations consistent with a flat- 
tened geometry the BLR, it must be differe nt than the one sug- 
geste d for Narrow Line Seyfert 1 (NLSvl: lOsterbrock & Poggg 
Il985l) and other Population A sources; for example if only for 
the simple reason that the minimum FWHM observed in face- 
on RL quasars is at about 3000 km s^^ broader than the nar- 
rowest Population A Balmer profiles ^ 500-600 km s~^, whose 
flattened BLR is also assumed observed face-on (e.g.. lZhou et al.l 
l2006l ; IZamfir. Sulentic & MarzianilbO OS*). The assumption of Kep- 
lerian motion in the BLR (and of FWHM(H/3) as a virial estima- 
tor) finds support from the discovery that the minimum observed 
FWHM(H/3) increases systematically with source luminosity (from 
^ 1000 km s"^ at logL^oKerg s"^) = 44 to ?a 3000 km s"^ at 
logL(,oi(erg s~^) = 48). This trend is easily explained if we are 
observing line emission from a Keplerian disk obeyi ng a Kaspi 
( lKaspietal.ll200dl2007h relation dMarziani et"ai]|2009l) . 

1.2 The Drivers of Type 1 AGN Diversity 

Numerous studies sought to explain the quasar spectroscopic di- 
versity in terms of physics (intrinsic to the AGN) coupled with 
orientation of the observer relative to the AGN structure and ge- 
omet ry. Much emp hasis was placed on properties like BH mass 
(e.g., lBorosonl2002|). luminosity (e.g. . Marziani et al. 200 9j) and ac- 
cretion rate (e. g., ^Sulenti c et al ] |2000b ; iKuraszkiewicz et al ] |2000l ; 
iMarzianiet 30 2001. 20033) as dominant drivers of the quasar di- 
versity. 

From a practical perspective the derivation of these quantities 
is complicated by several arguments, of which the most cumber- 
some are: 1) the nature of the underlying continuum is rather poorly 
understood, 2) numerous iron lines (Fell and Felll) contaminate the 
UV-optical-IR window producing also a pseudocontinuum, 3) not 
all emitted broad lines are signature s of virialized gas gravitation- 
ally bound to the central m as s (e.g., lGaskell|[l982l ; iMarziani et al.l 
1 19961 ; ISulentic et ai] |2007|) . 4) broad lines can show inflections 
and asymmetries th at may indicate a composite nature of the 
emitt i ng region (e.g.. IMarziani et ai] |l996l ; lNetzer & TrakhtenbrotI 
I2OO7I ; IMarziani et alj|2009h . likely connected to a complex struc- 
ture and kine matics of the emitting region even for the same ionic 
species (e .g., 'Sulen tic et al. I l2000bl ; iNetzer & TrakhtenbrotI |2007| ; 
lMarziani~ et al. 200^^, 5) physics and orientation are coupled and 
spectral information reflects both. 

Other contributors have also been suggested o r specu- 
lated ( mainly for RL/RQ): the B H spin ( Wilson & Colberll 
I1995I ; iModerski. Sikora & Lasotal 1 19981 ; IMeieil I2OOII ; 



iTchekhovskov. Naravan & McKirmev |200Si), the relatio n to 



the host galaxy (e.g., IWoo etal.1 12005; Oht a et alj 12007 ), the 
host galaxy morphology (e. g. Icapetti & Balamverdd 20061 ; 
ISikora. Stawarz & Lasota 2007) and/or its l ink with the nucleus 



(e.g., [Hamilton, Casertano & Tumshek' 200^ or the environment 
(e.g., Kauffmann, Heckman & Best 2008.). Some researchers 
prefer the orientation (either internal or ex ternal) for explainin g 
the differences between all type 1 quasars dRichards et"al]l2002l) . 
perh aps to remain more smo othly anchored to the unified model 
(e.g.. lUrrv & Padovani|[l995h . 

Metallicity was also fit into the picture of quasar spectroscopic 
diversity, s everal different correlations being reported: metallicity- 
luminosity dHamann & Ferl and'l999^;' Nagao, Marconi & Maiolind 
2006), metallicity-BH mass (Warner, Hamann & Dietrich 200^ 
and metallicity-accretion rate ( S hemmer et al.ll20040 . The strength 
of the optical Fell BLR emis sion may be connected with th e abun- 
dance/metallicity of quasars dNetzer & Trakhtenbrolll2007h . but all 
those results are difficult to relate and interpret in the framework of 
the Eigenvector 1 space for the following reasons: (1) BH mass is 
typically larger in Population B sources, which also show weaker 
Fell emission, (2) accretion rate is typically larger in Population A 
sources, which show more prominent/stronger Fell emission and 
(3) luminosity is part of Eigenvector 2, formally orthogonal to 
Eigenvector 1 space. The Eigenvector 1 h as also been interpreted 
as an evolutionary or "age" sequence (e.g., lGrupell2004 and refer- 
ences therein). Population A sources that are best described by low 
BH masses and high Eddington ratios may represent an early stage 
in the evolution of AGN, in that sense they may be young. 



1.3 The 4DE1-Based Approach 

The diversity of quasars has been explored from various directions: 
RL versus RQ, NLSy 1 versus BLSyl and Population A versus Pop- 
ulation B (e.g.. lBorosoj2002l ; lMarziani et al.l2003bl ; ISulentic et al.l 
2006). 4DE1 appears as a promising tool that organizes the spec- 
tral diversity of quas ars maximizing t he dispersion o f their prop- 
erties ( Sule ntic. Marziani & Dultzin-Ha cvan 2000a; Sulentic et al.l 
L2000b, 2QQ%. We explore quasar phenomenology in the context of 
the 4DE1 space which involves optical, UV and X-ray measures for 
low redshift sources. Pre-SDSS work was based on our own Atlas 



of quasars spectra, usually b righter than B~16.5 (Marzia ni et al 
19961; ISulenticT Marziani & Dultzin-Hacvan 2000a; Sulentic et al 



2000bl ; IMarziani et alj2003al f 



Empirical evidence (e.gT iSulentic et al]|2007l) argues against 
indiscriminate averaging of spectra. One requires some sort of con- 
text within which to define the averages such as the RQ-RL pop- 
ulations. As sample sizes and data quality improve it becomes 
possible to use ever more refined contexts. The optical plane of 
the 4DE1 uses two of the most fundamental parameters involving 
FWHM(H/3) and Kpeii- In the case of our bright SDSS DR5 sam- 
ple a composite average spectrum of all N=469 individual spectra 
(see next section) would give us little insight into the nature of the 
BLR region in quasars. Averaging spectra with FWHM(H/3) cov- 
ering a range from 1000 to 16000 km s~^ and RfeJ/ from to 
2.0 would be an oversimplification of the quasar phenomenon. The 
analogy to indiscriminate averaging of OBAFGKM stellar spectra 
is obvious. 

We have compared quasar spectra in the population A-B con- 
text where sources are separated at FWHM(H/3)~ 4000 km s~^ 



4 Zamfir et al. 



(Figure la). In the context of the present study we stress two im- 
portant facts: 

• NLSyl sources, customarily defined by FWHM(H/3) < 
2000 km s^^, do not form a distinct class of objects relative to the 
bulk of Population A sourceflfMarziani et al. 2001) if we consider 
their optical/UV sp ectra OVeron-Cettv. Veron & Goncalvesll200ll : 
ISulentic et aT]|2002h . 

• There is a remarkable change in the H/3 line profile at 
FWHM(H/3)si4000 km s"^ H/3 is fit by a Lorentzian function up 
to FWHM(H)3)«40 00 km s~^, while a double Gaussian is ne eded 
for broader sources jSulentic et al]|2002l : iMarziani et ai]|2009l) . 

It is interesting to note that if one uses the RQ-RL separation to 
define this boundary wit hin the optical plane of the 4DEI spa ce (see 
Table 3 and Figure 4 of IZamfir. Sulentic & MarzianilbOOSh then a 
2D Kolmogorov-Smimov (2D K-S) test yields FWHM(H/3)«3900 
km s^^. RQ-RL distinction appears less general and subordi- 
nated to the A/B framework since Population B so urces can be 
both RQ and RL feamfir. Sulentic & Marzianill2008h . We are not 
claiming that the limit at FWHM(H/3)«4000 km s"Ms abso- 
lute or fixed. Several lines o f evidence suggest a dichotomy at 
FWHM(H/3)«4000 km s'^ l lSulentic et al1f2007l : IMarziani et al] 
[2009J ; especially impressive is the H/3 profile change mentioned 
above (and also discussed later in the paper), although as yet we 
cannot exclude the possibility that there is a smooth change of prop- 
erties along the optical El sequence. 

The first result listed above suggests that Population A sources 
are simply an empirically-motivated extension of the NLSyls to a 
somewhat broader Umit. The limit at FWHM(H/3)?5i4000 km s"^ 
is also motivated by the rarity of sources with broader H/3 if Rfe// 
> 0.5. An underlying physical parameter could be the accretion 
ra te, and the limit may be dependent on luminosity (see Figure 1 1 
of IMarziani et alJl2009l for the locus of fixed accretion rate in the 
FWHM vs. Luminosity plane). In this case, the effect of orientation 
may blur the boundary at a critical accretion rate potentially asso- 
ciated with a BLR structure transition, creating an almost constant 
limit at 4000 km s^^ at low and moderate luminosity. No break 
is observed near FWHM(H/3)Ri4000 km s"Mn the 4DE1 optical 
plane and it would be truly remarkable if one were present given 
also measurement uncertainties. Summing up, several lines of ev- 
idence from previous work mandate the consideration of - at the 
very least - a limit at FWHM(H/3)~4000 km s~^, lest gross errors 
and misconceptions are introduced in the data analysis. 

The paper is organized as follows: § 2 presents a detailed char- 
acterization of H/3 broad emission line in low redshift SDSS spec- 
tra based on: composite median spectra in two representative bins 
defined within the optical plane of the 4DE1 parameter space and 
a set of diagnostic measures like asymmetry, kurtosis and centroid 
shift. It also explores the connection between these diagnostic mea- 
sures and the 4ED1 space. In § 3 we compare sources of low and 
high Eddington ratio both in terms of line measures and composite 
median spectra. § 4 and § 5 are dedicated to discussion and conclu- 
sions. Throughout this paper, we use Ho = 70 km s~^ Mpc~^, Q,m 
= 0.3 and Q.k = 0.7. 



^ A recent study iNikotaiuk. Czemv & Gurvnowicj )2009l) postulates that 
"the division between NLSyl and BLSyl based on the width of H/3 hne is 
less generic than according to the soft X-ray slope", which actually is one 
of the four adopted dimensions of the 4DE1 space. 



2 DETAILED CHARACTERIZATION OF H/3 EMISSION 
LINE PROFILES IN LOW Z QUASARS 

2.1 General Remarks 

The robustness of any systematic study of quasar broad emis- 
sion lines is conditioned by two things: a) availability of large, 
representative samples and b) the quality of spectra, both in 
terms of resolution and signal/noise (S/N) ratio. The situation 
is greatly improved thanks to dedicated surveys like the Sloan 
Digital Sky Survey (SDSS) that make available to the community 
large, homogeneous data-sets. This work involves exploita- 
tion of - 470 bright low z «0.7) spectra from SDSS (DR5; 
lAdelman-M cCarthv e t al] (20071). No comparable sample exists 
in terms of uniformity, completeness and spectroscopic quality. 
Construction and characteri s tics o f the sample are explained in 
IZamfir. Sulentic & MarzianH ( l2008h . Our work focuses chiefly on 
the two optical dimensions of the 4DE1 space, with occasional 
commentaries relative to the other two based on previous studies. 
One goal of this study is to emphasize that composite average 
spectra - averaged in a well defined context like 4DE1 space - can 
reveal the underlying stable components of the broad lines. We 
address or readdress a few fundamental questions in this paper: 

• Can we further emphasize quasar diversity employing measures 
such as centroid shifts, asymmetry and line kurtosis? 

• How are these diag nostics tied to and/or in corporated within the 
4DE1 space (see also lBoroson & Greeijl 19921) ? 

• Is the distinction between Population A/B more general than RL 
versus RQ quasars? 

• Does it reveal evidence for a Population A/B dichotomy or 
continuity across 4DE1? 

• Can we reveal/infer some kinematics/structural differences with 
the aid of the aforementioned line diagnostics? 

The departure point for our present analysis is the sam- 
ple of N=477 bright SDSS (DR5) quasars (z<0.7) defined in 
IZamfir. Sulentic & Marziani| j2008l) . The sample was defined such 
that psf g or i-band magnitudes are brighter than 17.5, representa- 
tive for the luminosity regime log[Li,oi(erg s~^)]=43-47. The red- 
shift limit was chosen so that we have a good coverage of the H/3 
and its adjacent spectral regions employed for the underlying con- 
tinuum definition. The brightness limit on the apparent magnitude 
was driven by the need of good quality spectra, which allows re- 
liable measures on the broad line profiles. Figures 1 and 2 present 
the general properties of the sample in terms of location in the opti- 
cal plane of the 4DE1 space, other spectral properties of the broad 
lines, the number of FRII and core-dominated RL sources along 
with the estimated radio loud fraction in each bin, and the S/N (sig- 
nal to noise) spectral quality, respectively. For a detailed presenta- 
tion of the criteria i nvolved in the construction of our sa mple please 
refer to section 2 of lZamfir. Sulentic & Marzianil ( I200S) . 

Figure la shows the source distribution in the optical plane 
of 4DE1 whose parameters involve broad H/3, FWHM(H/3), and 
the equivalent width ratio of the optical Fell (4570A blend) and 
broad H/3, RFei/=W(FeII4434-4684A)/W(H/3). Sources show the 
characteristic "banana" shape distribution found in earlier studies. 

2.2 Spectral Bins in the Optical Plane of the 4DEI 

Bins of AFWHM(H/3) = 4000 km s"^ and ARpeii = 0.5 
are shown in Figur e 1 following the definitions introduced in 
ISulentic et al] ( |2002|) . The width ARpeii = 0.5 basically corre- 
sponds to the measurement uncertainty at ±2a confidence level. 
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For example, if we consider tlie center of bin A2 (Rfb// = 0.75), 
data points are not significantly different from this central value 
at a confidence level ±2a since the error of individual measure- 
ments is in the range ±0.20 to ±0.25. A few extreme sources with 
FWHM(H/3) > 16000 km s~^ or Kpeii > 2.0 are not consid- 
ered here. The two sources that are located in what would be bin 
B3 are ignored at this time as well. The bulk of sources lie be- 
tween FWHM(H/3)= 1000- 12000 km s~^ and Rpeii < 1.5. Our 
sample does not include a significant population of NLSy 1 sources 
between FWHM(H/3)=500-1000 km s"\ which are classified as 
"Galaxy" in SPS S. Such SDSS sources have been tabulated by 
Zhou etal J (120061) and n=40 of them were included for evaluation in 
Marziani et alj l 2009). In our present report we make u se of the Fell 
optical template provided bv Boroson & Greenl ( ll992D study, which 
is not suitable for very narrow broad lines. In only 14 s ources (less 
than 3%) do we measure a negligible RfcII (see also toong et al.l 
[2009; Fell may be a ubiquitous signature in Type 1 AGN and we 
speculate that Fell presence may be incorporated in the empiri- 
cal definition of Type 1 phenomenology). The median centroid of 
the source distribution is located at (Rpeii, FWHM(H/3)) = (0.53, 
4300 km s~^). The distribution shows no obvious concentration at 
the centroid values, instead we see a broader concentration between 
O.I-l.O, 1000-6000 km s"\ The source distribution is "elongated" 
(like a banana) and largely reflects the intrinsic source dispersion 
rather than measurements errors (~ 20ct in FWHM(H/3) and ~ 16a 
in Rpeii)- 

A trend is seen in the sense that narrower FWHM(H/9) sources 
show a large dispersion in Rfe// and weaker Rpeii sources show 
a large dispersion in FWHM(H/3). The range ofRpeii is similar at 
FWHM(H/3)= 8000 or 12000 km s~^ and the range of FWHM(H/3) 
is the same at Rpeii = 1-0, 1.5 or 2.0. Sources with very broad 
FWHM(H/3) and large Rpeii measures are not present in this sam- 
ple of quasars and are likely either nonexistent or obscured. The 
distribution shows tails at both sides, i.e., extreme FWHM(H/3) (up 
to 40,000 km s"\ Wang et al. 2005) and extreme Rpeii (up to 
2.5); the tails are mutually exclusive and they represent about 5% 
of the total sample. 

Figures lb and Ic provide median and mean values for key 
parameters in each bin: the number of sources and average val- 
ues ofbolometri£jummos^^ lOXLx, where A = 5400A, 
see lZamfir. Sulentic & Marzianil2008l) . redshift z, equivalent width 
of H/3 and of Fell (4570A blend), FWHM(H/3) and Rpeii- These 
numbers are representative for quasars of intermediate (bolomet- 
ric) luminosity (logLboi[erg s^^] = 43.5-46.5) at low redshift (z < 
0.7). RL sources tend to strongly prefer the Population B regime. 
Nonetheless, there is a small fraction of RL amongst population A 
sources, but we should also emphasize that ~40% of them (7 out 
of 18) show FWHM(H/3) in very close to the 4000 km s"^ A/B 
adopted boundary. 

Comparing the average (median or mean) equivalent width 
values W(H/3) and W(FeII)(Figures Ibc) we see an increase of 
about 45% in the former and a decrease of about 50% in the lat- 
ter measure in bin Bl relative to bin A2. Rpeii consequently de- 
creases by a factor 2.8 from A2 to B 1 . RL quasars strongly concen- 
trate above FW HM(H/3) ^ 4000 km s~^ and belou' Rpeii ~ 0.50 
(see Figure 4 of IZamfir. Sulentic & Marzianill2008l) . Most excep- 
tions (i.e.. Population A RL) show core dominated radio morphol- 
ogy suggesting that they are viewed near the radio jet axis. If H/3 
emission arises from an accretion disk oriented perpendicular to the 
jets' axis then a correction to FWHM(H/3) would move most core- 
dominated RL from Po pulation A into the region of Population B 
jMarziani et"ai]|2003bl) . RL sources showing lobe-dominated FRII 



and core-dominated radio morphologies sho w a clear and signif- 
icant separation in the 4DE1 opt ical plane i lSulentic et a"i]|2003l ; 
IZamfir. Sulentic & Marzianill200g ). This is a hint that the orienta- 
tion plays an important role at least for RL quasars and possibly for 
all quasars. 

A more extreme comparison (relative to Figures Ib-c) could 
involve the less populated bins A3, B2 and B1+, each with about 
30-40 sources. Average luminosities are very similar. On average 
W(Fell) in bin A3 is about 4x larger than in bin B1+. Average 
W(H/3) in bin B1+ is about 1.5x larger than in bin A3. There are 
no RL sources in bin A3, but a large concentration of RL is found in 
B1+ (see Figure 4 in lZamfir. Sulentic & Marzian i 2008). We avoid 
bin Al with the smallest FWHM(H/3) and Rpeii measures because 
it may involve a mixture of Population A and B sources. In terms 
of spectral quality we estimate mean and median ratio S/N ~ 29 
(defined in the continuum region 5600-5800 A). Fi sure 2 shows 
the (histogram) distribution of pixel S/N ratios for the N=469 SDSS 
spectra in our sample. 



2.3 Composite Spectra of Quasars 

Composite average quasar spectra have proven to be a pow- 
erful tool for revealing both commonalities among and dif- 
ferences between quasars. High S/N spectral composites have 
served many purposes. First they allowed a general descrip- 
tion of the quasar phenomenon by: a) providing a refer- 
ence or cross-correlation template spectrum by identifying (and 
measuring) weak emission features (strengths, widths, shifts, 
asymmetries, etc.), otherwise hardly observed in individual 
spect ra (e.g., IWills etal. 1980; Veron-Cetty, Tarenghi & Vero^ 
| l983l; Ibo^ Il990l; iFrancis et alj 119911; lzhengetal.1 Il997l; 



Ivanden Berk et al.l l 200ll; IConstantin & Shieldsl |2003| ; IScott etaU 
2004; Glikman. He lfand & White! |2006|) ; b) defining the un- 
derlying continuu m by is olating wavelength regi ons (almost) 
free of lines (e.g., Francis et a l. 1991; Vanden Berk et al. l 1200 ll ; 
IConstantin & Shieldsll200'3l) r c) revealing the omnipresence of Fell 
bl ends throughout the UV an d optical domains (see e.g., section 4.2 
in IVanden Berk et alj|200ll and references therein); d) evaluating 
colors of quasars and/or K(z) corrections jCristiani& ^119901 ; 
Iblikman, Helfand & White 2006). 

As larger samples of quasars spectra became available, com- 
posites were constructed for sources with very similar properties 
like: 1 ) radio-loudness (e.g., [Cristiani & Vio 1990; Zheng et al 
1997; Brotherton et al. 2001; Telferetal. 2002; Bachev et al 
,2004; Glikman. Helfand & White 2006; Metcalft & Magl iocchetTi 



l2006l ; Ide Vries. Becker & Whitell2006h. 2) radio core-to-lobe flux 
density ratio (e.g.. Baker & Hunsteaj 1995|). 3) radio morphologv. 
FRII versus core-dominated, e.g.^ Ide Vries. Beck er & Whi2 
(2006), 4) optical continuum slope or relative color jRTchards et aP 
2003l). 5) broad emission line velocity shifts teichards et alj|2002l: 
Huetal.ll2008bh. 6) CIV A1549A line width dBrotherton et al. 



1994[) . 7) emission line wi d ths a nd strengths (Sulentic et al 
2OO2I ; IConstantin & Shieldsl |2003| ; iBachev et al.l |2004|) . 8) 
Black Hole Mass ( Mbh) and Eddington ratio (Lj,oi / LEdd) 
(e.g.. iMarziani et al.l l2003bh . 9) variability dWilhite et tJI |2005| ; 



Perevra et al. 20061). 1 0) Broad Absorptio n Lines - BAL (e.g 



Brotherton et alj|200ll ; iRichards et aL )'2002'; 'Reichard et al.l l2003h 



and II) soft X-ray brightness (e.g.,[Green 1998). Such compos- 
ites facilitate an emphasis on differences between properties of 
quasars. The luminosity and redshift dependence/evolution of 
quasar prope rties was sought and tested with the aid o f composite 
spectra (e.g., ICristiani & ^Il990l ; IZheng etalj|l997l ; iBlair et al.l 
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20001; Croom et al. 20021: ICorbett et alj 20031; Yio et alj 2004 


Naeao. Marconi & Maiolind 20061; Scott et alj 2004; 


Bachev et al.l 


20041; iMetcalft & Ma£liocchetti"2006l;lMarziani et al. 


l2009h. 



An important que.stion involves how to combine individual 
spectra into a composite. Moreover, how well does a composite 
spectrum preserve the properties of the input components and to 
what extent is it representative of the sample from which it has 
been constructed? It has been noted that there is no unique best an- 
swer to such questions. While a median composite preserves the 
relative fluxes of the emission features, a geometric mean compos- 
ite preserves better the shape of the underlying continuum (e.g., 
IVanden Berk et a i] |200ll) . Additional questions can be related to 
how to weight each spectrum's contribution into the composite. De- 
pending on the number and quality of individual spectra available 
one could decide whether a weighte d -average (arithmetic mean) 
composite (e.g., lFranci setal]|l99ll ; iGlikman. Helfand & Whit3 
[2OO6) is more usef ul than a composite where a ll spectra are given 
equal weig ht (e.g.. lConstantin & Shield"sll2003[) . The most impor- 
tant thing is the context within which we define the composites; we 
consider 4DE1 is the best approach available. 

Given the S/N distribution (Figure 2) of our sample, the homo- 
geneity of the SDSS spectra and the nature of our spectral analysis 
we see the median composites with equal weight from individual 
components as most suitable (for a disc ussion on the releva nce of 
the median composites see section 4.1 of lBachev et al.|[2004h . 



2.4 H/3 Profile Shape 

2.4.1 Data Analysis 

Low z rest-frame spectra require galaxy decontamination and flux 
corrections. They are normalized to the 5100 A specific flux prior 
to being considered for the composite. The full analysis procedure 
for a quasar optical spectrum focusing on the H/3 region is outlined 
in|Zamfir, Sulentic & Marziani (2008). We added an intermediate 
step in the present study by defining the underlying continuum as a 
simple power-law anchored in the regions 4195-4218 A and 5600- 
5800 AThe Fell te mplate used to model Fell emission around H /3 
is based on IZwl (Boros on&Greedll992l ;l Marziani et al ] |2003al) . 
The template matches the 4750 A Fell "blue" blend significantly 
better than it does the Fell blend redward of [OIII]A5008A ("red" 
blend). In many cases the "red" blend of Fell may be overestimated. 
In order to alleviate this problem we opted for an additional itera- 
tion in the reduction process of the composite spectra: 1) we define 
the underlying power-law continuum, 2) subtract the continuum, 
3) find the best match between the Fell template and the spectrum 
(adjusting the width and intensity of the template), focusing on the 
best match for the blue blend of Fell, 4) subtract this template so- 
lution from the original rest-frame spectrum (available before step 
1) and 5) redefine a power-law continuum using the wavelength 
regions 4195-4218 A and 5100-5300A. This continuum is slightly 
steeper than the initial choice. We now repeat steps 2 and 3. This it- 
eration significantly improves the match between template and the 
red blend of Fell. 

From the spectrum thus obtained (continuum and Fell-free) 
we cut out the wavelength range bracketed by HeIIA4687A and 
[OIII]A5008AWe simultaneously fit all lines in the selected inter- 
val, allowing for narrow components (NC) in the case of Hell, H/3, 
[OIII]A4960A and 5008 A lines. The [OIII] lines require an addi- 
tional component, in most cases blueshifted relative to the quasar's 
restframe indicated by the NC. We do not require a match in width 
between the NC of H/3 and the core (the narrower of the two) com- 



ponents of the [OIII] lines. We constrain the [OIII] 5008/4960 flux 
ratio R=3 allowing for a maximum 15% error. All components are 
modeled with Gaussian functions in the first stage. Subsequently, 
we test whether a Lorentzian function provides an improved fit to 
the H/3 broad component (BC) in terms of reduced and by visual 
inspection of residuals. All other emission lines are modeled with 
Gaussian functions. We use the Levenberg-Marquardt algorithrr0 
to adjust the parameter values in the iterative procedure. 

2.4.2 Composite Spectra in 4DE1 

Figures 3 and 4 present composite spectra for bins A2 and Bl, re- 
spectively. We show both the region (upper) from 3500-7000A with 
continuum+line model and an enlarged view (lower) of the H/3 re- 
gion from 4600-5 lOOA with details of the line model. We focus 
on bins A2 and B 1 as representative of source differences between 
Populations A and B. We avoid bin Al because it likely contains a 
mix of both Population A and B as well as because it shows a range 
of Kpeii values that is too similar to bin B 1 sources. 

Extreme A bin composites (A3 and A4) tend to have stronger 
Fell emission while extreme B bins (B1+, B1++) involve sources 
with broader H/3 profiles and weaker Fell emission. The systematic 
weakening of Fell between bins BI and B1++ (see also Figures 
lb,c) may be caused by the increase in FWHM for Fell lines that 
increases with FWHM(H/3). As the Fell blends broaden they create 
a pseudocontinuum that may lead to systematic underestimation of 
Fell strength. However, thanks to the high S/N of our SDSS sample. 
Fell is detected for the first time in most sources with FWHM(H/3) 
> 8000 km s~^, revealing that this effect is not as influential as 
feared. 

Bins A2 and Bl are well matched in terms of number of 
sources as well as mean/median redshift and logLbo! (Figures 1- 
b,c). N=129 and N=131 sources were used to generate the A2 and 
B 1 composites respectively. The differences between the two com- 
posites would be slightly larger if we included in A2 t he NLSyls 
source s with FWHM(H^ < 1000 km s"^ tabulated bv lZhouetall 
( I2OO6I) . However inclusion or exclusion of these sources will make 
no significant change in the results. Marziani et al. (2003b) show 
that sources with FWHM(H/3) < 2000 km s"\ often described 
as a separate class, are similar in terms of line profile shape (i.e., 
Lorentzian) to those having FWHM(H/3) in the range 2000-4000 
km s^^. 

The A2 median composite shows high enough S/N to reveal a 
modest inflection between the BC and NC of H/3. The NC adopted 
by the model shows FWHM ^ 350 km s~^, very similar to the 
core components of the [OIII] lines. BC shows FWHM(H/3)~2600 
km s^^. The best fit (Figure 3b) to the BC of H/3 profile is ob- 
tained with a Lorentzian function, confirming results o f previous 
studies jSulentic et al.l2002l ; lMarziani et alj2003bl.l2009l) . The pro- 
files of NLSyl sources have long been described as Lorentzian 
dVeron-Cettv. Veron & Goncalvesll200ll) . but that profile descrip- 
tion now extends to FWHM(H/3)«4000 km s"\ A Gaussian fit 
does not provide satisfactory results in light of residuals. A 
pair of unshifted Gaus sians provides a reasonable fit dPopovic et al.l 
l2004l ; lBon et al.|[2009l) since they essentially represent an approxi- 
mation of a Lorentzian. A Lorentzian is preferred as a description 
of the A2 H/3 profile based on x'^- The broad Lorentzian component 
shows no shift relative to the quasar's restframe. 

There is a positive residual redward of [OIII]A5008A. This 

3 www.nr.com: |Press, Flannerv & Teukolskvl fl986l) . 
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residual can be an artifact produced by a mismatch between tlie Fell 
template and the composite spectrum coupled with an oversimpli- 
fied definition of the underlying continuum. The red-shelf has been 
pointed out in Marziani et al. (2003a) (see their section 4.3) and 
IVeron. Goncalves & Veron-Cettvl ( |2002|) . The latter study focused 
on two sources with such a prominent feature and showed that the 
main contributors to the "shelf" are most probably the Hel lines AA 
4922, 5017 A and is not caused by a broad redshifted component of 
H/3 or [OIII] lines. Even if a VBLR is assumed in our case here and 
affects or causes the red "shelf", its contribution to the total flux of 
H/3 is no more than 10%, which is within the uncertainty range of 
the fitting routine. 

The median composite for bin Bl (Figure 4) shows striking 
differences from the A2 composite, most notably in the shape of 
the broad H/3 profile. We see two inflections: 1) between the NC 
and the BC and 2) on the red side of the BC. The simplest fit to the 
broad component therefore requires a double Gaussian involving 
a broad (BC) unshifted and a very broad (VBC) component, red- 
shifted by ~ 1500 km s~^, which also involves about 53% of the 
total H/3 broad line flux. This model is again in agreement with our 
earlier studies. Detection and measurement of the VBC requires 
spectra with good S/N and moderate resolution - the median com- 
posite suggests that it is very common if not ubiquitous in sources 
with FWHM(H/3)> 4000 km s^^. It is obvious that some or many 
sources show FWHM(H/3)> 4000 km s~^ because of this second 
component. If we interpret the unshifted BC as the classical BLR 
then we might expect it to show a Lorentzian shape. Fitting the B 1 
composite with an unshifted Lorentzian plus redshifted Gaussian 
yields a poorer and a VBC comprising less than ~ 10% of the 
total broad line flux, with ~ 6000 km s^^ redshift. Past and present 
results suggest that a double Gaussian is the best choice. The clas- 
sical "unshifted" BC Gaussian shows a only small blueshift of ~ 
-100 km s"^ and FWHM of ~ 4000 km s"\ 

The width of the BC component in the bin Bl composite is 
significantly less than the median value of FWHM(H/3) ~ 5400 
km s~^ reported for bin B 1 in Figure lb. We conclude that FWHM 
measures of broad H/3 that do not correct for the VBC will lead 
to serious BH (Black Hole) mass overestimates. At the same time, 
for the bin A2 there is no difference between the median value re- 
ported in Figure lb and the measured FWHM(H/3) in the compos- 
ite median spectrum. This suggests that no VBC correction is re- 
quired for Population A sources making their H/3-based BH mass 
estimates more reliable. The bin A2-B 1 comparison suggests fun- 
damental structural and kinematic differences between Population 
A and B sources (see also Sulentic et al. 200Z) most likely coupled 
with non-trivial effects of orientation. 



single source. Since many of the features do not show up in the 
composite they presumably represent relatively short-lived, tran- 
sient and therefore lees important changes in the profiles. 

We want to go beyond the optical plane of the 4DE1 space 
because we find that sources inside bins like A2 and Bl, while 
showing similar FWHM(H/3) and RfeI/, also show considerable 
profile diversity implying that the two measures, width and relative 
strength of Fell emission, are not enough to characterize the diver- 
sity. For each spectrum shown in Figure 5 we also indicate the bin 
in the optical plane of the 4DE1 space, the bolometric luminosity, 
spectral continuum S/N (as defined earlier) along with measures of 
centroid shifts, asymmetry and line shape (see next section). This 
set of 32 spectra is assumed representative in terms of S/N ratio, 
luminosity and occupation of the 4DE1 optical plane. In addition 
to the already stated purpose of presenting the great variance of 
the Balmer H/3 line properties amongst our low z/intermeidate lu- 
minosity sample, this subset serves also the purpose of estimating 
typical uncertainties for the line diagnostic measures described in 
detail in the next section. 



2.4.4 Higher Order Moments of HP Profiles as Diagnostics of 
the BLR Physics 

In order to get a deeper insight into the hypot hesized ex istence 
of two quasar popu lations A and B (e.g., Sulentic et al. 20Q3; 
iMarziani et alj200 9') we examine the sample distributions of higher 
order moments of the H/3 profiles. We focus on centroid shifts, 
asymmetry and kurtosis (line shape). The H/3 profile is fit with 
a high order polynomial making no assumptions about different 
components that might be present. Higher order moments of the 
broad lines are affected by larger uncertainties than the line width 
measures, nonetheless, they offer the possibility to further constrain 
physical models of the BLR. We benefit in this analysis from the 
uniformity of SDSS spectra in S/N and resolution. 

For quantitative definitions of the parameters involved i n our 
study we consider those proposed in IMarziani et alj ( Il996l) . For 
clarity we reproduce them here in the following order: Asymme- 
try Index A.I., Kurtosis Index K.I. and centroid shift c{x): 



A.I. 



Ur,fl(l/4) + Wr,B(l/4) - 2Vr,peak 



K.I. 



^Jr,fl(l/4) -w,,s(l/4) 
«r,fl(l/4) -u,,s(l/4) 



(1) 



(2) 



2.4.3 Diversity of H(5 Profiles of Individual Spectra 

We have argued that indiscriminate averaging of spectra will ob- 
scure the most physically fundamental differences in a source sam- 
ple. The purpose of this section is to go to the other extreme and ex- 
amine source spectra individually. Figure 5 shows a selection of in- 
divid ual H/3 pro files in order to illustrate that all previously known 
(e.g., IStirpel 199Q) diversity is present in the SDSS sample. The A2 
and B 1 median comp osite spectra present ed and discussed earlier 
(and also discussed in ISulentic et aLl l2002') do not show the variety 
of features observed in individual, single-epoch spectra shown in 
Figure 5. This justifies the interpretation that average spectra rep- 
resent the underlying commonality, more stable components in the 
H/3 profiles of different sources. The individual spectra are, to some 
extent, equivalent to the results of spectroscopic monitoring of a 



The \r,R and v,.,b refer to the velocity shift on the red and 
blue wing, respectively (at the specified hight of the profile), cal- 
culated relative to the restframe of the quasar. By "peak" we mean 
9/10 fractional intensity. The centroid shift will refer to x = 1/4 
fractional intensity (what we call the line "base") or a; = 9/10 (the 
line "peak"). We point out that all these measures avoid the region 
below 1/4 fractional intensity due to the larger uncertainties asso- 
ciated with the definition of line wings and underlying continuum 
plus Fell pseudo-continuum. 

2.4.5 Higher Order Moments and the Optical Plane of the 4DE1 
Parameter Space 

It is important to identify any kind of relation between 4DE1 
and higher order diagnostic parameters that characterize the pro- 
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file of the H/3 profile. Figure 6 shows plots of FWHM(H/3) versus 
the newly defined diagnostic measures: (a) vs. asymmetry index 
A. I. (1/4), (b) vs. kurtosis line shape measure K.I., (c) vs. normal- 
ized base shift c(l/4)/FW025M and (d) vs. normalized peak shift 
c(9/10)/FW0.90M). Whenever shown, the two vertical lines indi- 
cate the 2o uncertainty lever on either side of zero. We also indi- 
cate the two populations A and B defined in terms of FWHM(H/3) 
less than and larger than 4000 km s~^, respectively. In each panel 
we also display the FRII radio sources with distinct symbols (blue 
solid). A few commentaries may be worthwhile: 1) There seems to 
be a trend in panel (a) in the sense that lower widths profiles are 
more symmetric, but sources with broader profiles show a numer- 
ical excess of (larger than the 2o uncertainty) of red-asymmetric 
basis; 2) Panel (b) shows no clear trend, but it indicates that the 
kurtosis values of broader H/3 profiles cover more uniformly a 
wider range of vales than the narrower profiles; 3) Panel (c) re- 
veals bears some similarity to panel (a), with a subtle difference; at 
lower widths there seems to be a slight excess of blueshifted "ba- 
sis", and the slight excess of redshifted basis is still apparent for 
broader profiles; 4) there is clear excess of sources with blueshifted 
"peaks" over the whole range of FWHM. FRII sources do not seem 
to have any special behavior when contrasted against the other Pop- 
ulation B sources. 

The centroid shift measures are normalized in each case to 
the width of the profile at that same fractional intensity. We dis- 
tinguish in each case Population A and B sources, as defined in 
the context of the 4DE1 parameter space. The vertical dotted lines 
indicate the 2a median uncertainty derived for the subset of 32 
sources shown in Figure 5. Individual source profile uncertainties 
(for the 32 spectra) are estimated allowing for ± 5% variation in 
fractional intensity and then propagating the errors into the final di- 
agnostic measure. Thus, the typical errors are: 2a(A.I .[1/ 4)) ~ 
0.16, 2a{c{l/A)/FWQ.2hM{Hp)) ~ 0.07, 2a{K.I.) ~ 0.092, 
2a {c{9/ 10) /FW0.90M{H 13)) ~ 0.22, 2a (c(l/4)) ~ 520 km 
s~^ and 2a (c(9/10)) ~ 320 km s"\ 

Figures 7a-d show the histogram distributions of the asym- 
metry index A.I.(l/4), the centroid shift at 1/4 and 9/10 fractional 
intensities c(l/4) and c(9/10), respectively and the kurtosis measure 
K.I. The width of the bins in Figure 7 (and 8) are chosen such that 
they are approximately equal to la typical (i.e., average) uncer- 
tainty. In Table 1 we present the average (mean and median) values 
of the parameters shown in Figure 7. We make the following com- 
ments relative to Figure 7: i) Population B sources show an excess 
(outside 2a) of red asymmetric profiles and profiles red-shifted at 
the base, ii) Population A sources are more symmetric and they 
show a slight excess of blue-shifted profiles at 1/4 from zero-level, 
iii) there is an excess of blue-shifted peaks in both Population A 
and B sources and iv) Population B sources show an extended tail 
of large kurtosis values, while Population A quasars show a more 
symmetric distribution about K.I.=0.35, with a clear preference for 
lower K.I. values. A one dimensional K-S test indicates a 4% prob- 
ability of null hypothesis when comparing the normalized "peak" 
shift distributions (lower left panel of Figure 7). For all other mea- 
sures shown in Figure 7 a ID K-S test shows a null-hypothesis 
probability less than 10~*. The numbers in Table 1 appear to con- 
firm the aforementioned observations i) and ii). 

Figure 8 shows a comparison of "base" and "peak" shift for 
Populations A and B, this time without normalization to the line 
width at the corresponding fractional intensity as in Figure 7. We 
are aware that the un-normalized centroid shift values are less 
meaningful than the normalized ones. Nonetheless, the two pan- 
els in Figure 8 seem to show: 1) the two populations are signif- 



icantly different (Population A sources being distributed mostly 
within ±2cr, while Population B quasars showing large excesses of 
blueshifted peaks and redshifted basis outside the 2a range and 2) 
the measured shifts are quite large, "base" redshifts over 1500-2000 
km s^^ and "peak" blueshifts in excess of -1000 km s^^ being ob- 
served. 

We tested whether the observed distributions are spurious ef- 
fects of the S/N. For this task we reconstructed the histogram dis- 
tributions shown in Figure 7 restricting the sample to sources of 
S/N < 15 20% of our whole sample; see Figure 2). All dis- 
tributions presented above are confirmed and are not caused by 
the spectral S/N. We also tested if the distributions presented in 
Figure 7 are driven by the preferential addition of a subset of 
SDSS quasars selected with some restric tion on their radio prop- 
erties jZamfir. Sulentic & Marzianill2008l) . For this reason, we iso- 
lated only the RQ sources of our sample (i.e., N=388 sources with 
logLa Icm < 31.6 erg s"^ Hz" \). Although we don't show the new 
histograms here, the distributions are qualitatively very similar to 
those in Figure 7, thus we have clear indication that the population 
A/B concept is more general t han the RL/RQ dist i nction . 

The pioneering study of iBoroson & GreenI ( Il992h incorpo- 
rated measures of asymmetry, line shape and shift in their set of 
quasar spectral measures. They reported a significant correlation 
between Rpeii and the asymmetry of H/3 (see their Figure 5). 
Sources with strong Fell emission tend to show a flux excess on the 
blue wing of H/3 and the weaker Fell emitters tend to be show more 
red-asymmetric H/3 profiles. The larger sample under investigation 
here reenforces their result (Figures 9ab), the trend being clearly il- 
lustrated in diagrams that combine A.I.(l/4), c(l/4)/FW0.25M and 

Rpell- 

Figure 9 presents the most interesting plots involving 
pairs of measures from the following set: Rpeii, A.I. (1/4), 
c(l/4)/FW0.25M, K.I., Lboi/^Edd- The plots involving the peak 
shift are not shown because they do not illustrate any particularly 
special aspects beyond the already mentioned ones (related to the 
histograms in Figure 7). 

Figures 9ab (log-linear plots) show the correlations between 
4DE1 parameter Rpeii and the asymmetry and "base" line shift 
measures. We also identify sources showing FRII radio morphol- 
ogy symbols (dark blue). The FRII quasars also serve as an inde- 
pendently (radio) defined subsample of Population B sources. Both 
panels show a trend from blue asymmetries/shifts at large Rfe/i 
values to red asymmetries/shifts at lower Rpeii values. A Pearson 
correlation test applied show to the whole sample yields a coef- 
ficient of rp = - 0.50 (with a probability of not being correlated 
P ~ 0) for panel (a) and r„ = - 0.58 (P ^ 0) for panel (b). The 
bisector fit ( Ilsobeetal.ll990h for (a) corresponds to the linear rela- 
tion A.I.(l/4)= -0.541ogRFe//-0.15 and for (b) c(l/4)/FW0.25M = 
-0.271ogRFeji-0.09. 

It is interesting to point out that in both cases the linear regres- 
sion indicates that when the c(l/4) or A. I. (1/4) yield a null value, 
the corresponding Rpeii ~ 0.5. The majority of FRII sources lie 
below the bisector fits and show positive A.I. (1/4) and c(l/4) val- 
ues. This indicates that Population A sources are less asymmetric 
and less redshifted (at profile "base") than Population B quasars 
(see also Figure 20 and related discussion). A 2D K-S test for RL 
vs. non-RL quasars within the optical plane of the 4DE1 space 
(Zamfir, Sulentic & Marziani 2008) indicates the same Rpeii '~ 
0.50 as a relevant boundary. 

Figure 9c shows the plot of Rpeii versus the line shape K.I. 
measure. We see a complex distribution, quite different from Fig- 
ures 9ab. We see different trends on either side of RFeij'~0.5. 
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Table 1. Average values for line asymmetry, kurtosis, base and peak shift for the two populations A and B, relative to Figure 7. 



Parameter 


Population A (N=260) 
mean (SD) ; median 


Population B (N=209) 
mean (SD) ; median 


A.I.(l/4) 


0.002(0.118) ; 0.001 


0.096 (0.155) ; 0.083 


K.I. 


0.356 (0.063) ; 0.350 


0.383 (0.092) ; 0.374 


c(l/4)/FW0.25M 


-0.011 (0.062) ; -0.010 


0.029 (0.076) ; 0.024 


c(9/l())/FW0.9()M 


-0.064 (0.277); -0.042 


-0.102 (0.300) ; -0.080 



The parent population of RL quasars (FRII sources) seems to 
prefer almost exclusively the region to the left of Rpeii = 0.50. 
Thus, we may already have several lines of evidence that seem 
to suggest RfcJ/ = 0.50 as a viable alternative to the original 
FWHM(H/3) = 4000 km s " ^ Population A-B boundary 

Previous quasar studies dSulentic et alj l2000bl : 
iKuraszkiewicz et al.l 1200 0*: 'Marzi ani et alj 1200 ll, l2003b l) indi- 
cated that the principal physical driver of the variance in 4DE1 
was the Eddington ratio htoifLEdd- In Figure 9d we show the 
correlation between logRpeii and logiLtoi/LEddY, l-Edd ~ 1-3 x 
10^** (Mbh/Mq) erg s^^ and 



Mbh{Mq) PS 5.48x10^ 



ALa(5100^) 



FWHM(HI3bc 



WOOkms- 



_ 1044 ergs - 

(see lSulentic et alj20o3) . 

A Pearson test yields r^ = 0.56 (P ~ 0) and a Spearman cor- 
relation test yields Ts = 0.59 (P ~ 0). We show the bisector fit 
(correlation coefficient p = 0.58). The linear fit logihtoifLEdd) 
= 1.471ogRF(;j/-0.32 indicates that Kpeii = 0.50 corresponds to 
an Eddington ratio htoi/l^Edd = 0.18. Fell emission increases 
in strength relative to H/3 as the Eddington ratio increases, i.e., 
higher accretion rate. We are now in a position to test whether 
the base asymmetry and centroid shift measures A.I.(l/4) and 
c(l/4)/FW0.25M are correlated with the accretion rate measured 
by hboifLEdd- Figures 9e-f show the two line diagnostic measures 
versus \og(hi,oifLEdd)- For the plot involving the A.I. (1/4) (panel 
e) a Pearson test gives rp= - 0.39 and a Spearman test gives rs = - 
0.43 (P ~ in each test). For panel (f) ip = - 0.28 (P ~ 7 x 10"^'^) 
and Ts= - 0.32 (P ~ 6 x 10""). The FRII sources prefer the low 
accretion regime. It appears that both the A.I. and the normalized 
c(l/4) measures are intimately connected to the 4DE1 parameter 
space. Nonetheless, their inclusion in this space definition may not 
necessarily lead to an extension of the number of dimensions. One 
of the two indices may be regarded as a meaningful surrogate for 
FWHM parameter for example, although the problem requires fur- 
ther investigation. 

Figure 9(g) shows the plot of the Eddington ratio vs. the line 
shape measure K.I.; the large majority of FRII sources are low 
accreting sources and their kurtosis measures spread over a wide 
range of values. It also seems to indicate that high Eddington ratios 
and and large kurtosis values are mutually exclusive. 



3 LOW VERSUS HIGH Lbol/Ledd SOURCES 

Several recent art i cles iHu beny et al. 2000; Marziani et al. '2003b'; 
CoUiri et aljr2006l;lBonning et al.i,2007. ; .Kelly et al...200& ; ,Hu et al., 
2008al ; iMarconi et al J 120091) came from various directions to the 
conclusion that around some critical hboifLEdd ~ 0.2±0. 1 the 



properties/structure of the active nucleus may change fundamen- 
tally, thus providing support for the notion that two distinct pop- 
ulations of quasars exist (Collin et al. 2006; Sulentic et al. 2007). 
While the concept of two Populations A and B is deeply anchored 
int o a long series of em pirical studies (see section 3 and Table 5 
of lSulenticetalJ|2007l for a summary), the nominal boundary at 
4000 km s~^ FWHM(H/3) may not have a very intuitive or imme- 
diate physical meaning or may even appear arbitrary. However, it is 
solely intended as a critical empirical measure that could indicate 
a fundamental change in the structure, geometry and kinematics of 
the BLR. 

n 2 It has been suggested (as was emphasized earlier) that the Ed- 
dj^il^ton ratio 'Lboifi^Edd is the main physical driver of the Eigen- 
vector 1 parameter space variance. Our sample reenforces these ar- 
guments (e.g., Figure 9d). Figure 10 compares the histogram distri- 
butions of bolometric luminosity L(,oi, black hole mass Mbh and 
Eddington ratio 'Lboif^Edd for populations A and B. Taken at face 
value they suggest that the two populations have very different dis- 
tributions of Msir and Eddington ratio, even though they show 
similar luminosity distributions (see Table 2). Obviously, given that 
the two populations have similar luminosity distributions, the sta- 
tistical difference in BH mass and Eddington ratio is a natural and 
direct consequence of the line width alone. Populations A and B 
being defined based on the FWHM(H/3). A ID K-S test in each 
case indicates a 1.2% null-hypothesis probability for the bolomet- 
ric luminosity distribution of the two populations and a probability 
less than 10~* for the other two parameters. However, although the 
AGN central mass and the Eddington ratio may be more physically 
meaningful quantities, FWHM is a more direct and robust measure, 
much less affected by uncertainties and derived with fewer assump- 
tions. 

We decided to construct two test samples which would 
be representative of high and low accretors; all sources with 
Lboi/LEdd > 0.3 (N=123) and with Lboi/LEdd < 0.1 (N=177) 
are assumed to be high and low accretion quasars, respectively. This 
distinction turns out to be representative of the previously adopted 
populations A and B. All high hboifl^Edd in our definition are part 
of population A and 87% of the low hboifLEdd are population B, 
only 13% spilling over the 4000 km s^^ boundary into population 
A. Obviously we are aware that this kind of distinction would auto- 
matically introduce some bias, i.e., the low accretors would prefer 
larger BH masses and less luminous sources and the high accretors 
would populate more the high luminosity and low BH mass regime. 

We point out a few remarkable results in Figures 1 la-d. Panel 
(a) shows source luminosity as a function of the normalized cen- 
troid shift at line "base" (1/4 fractional intensity). The bias indi- 
cated previously is clearly apparent. We may however focus on the 
region logLf,oj[erg s~^] > 45.5, where both high and low accre- 
tor sources are well represented. There is a clear separation of the 



10 Zamfir et al. 



Table 2. Average values for bolometric luminosity, BH mass and Eddington ratio for the two populations A and B, relative to Figure 10. 



Parameter 


Population A (N=260) 


Population B (N=209) 




mean (SD) ; median 


mean (SD) ; median 


log[L6oi(erg s"^)] 


45.4 (0.9) ; 45.6 


45.6 (0.8) ; 45.8 


log[MBif(M0)] 


7.8 (0.7) ; 7.9 


8.7 (0.6); 8.8 


^boll^Edd 


0.36 (0.28) ; 0.28 


0.07 (0.05) ; 0.06 



two samples in terms of base shift, with high accretors showing a 
preference for blueshifts and the low accretors for redshifts. Inside 
the ±2ct uncertainty interval there is some overlap while outside 
(most significant "base" line shifts) the separation is almost total. 
A few FRII sources show blueshifted profiles at 1/4 base and very 
few high accretors spill over the +0.07 boundary into the redshifted 
region. 

Panel (b) shows the normalized "peak" shift (at 9/10 fractional 
intensity) versus the asymmetry index A.I.(l/4). The separation is 
again quite striking: high accretors tend to show unshifted peaks 
with an excess of blue-asymmetric bases, while the low accre- 
tors distribute more toward blueshifted peaks and red-asymmetric 
bases. The simplest accretion d isk models predict redshifted profile 
bases and blueshift ed peaks dChen, Halpern & Filippeiikol Il989l : 
ISulentic etalll 19901 but see also lLaoij|l99l[ for X-ray line model 
predictions). We may find here support for the idea that low ac- 
creting sources produce their broad emission lines (the LIL/Balmer 
components) in a disk-like geometry. A 2D Kolmogorov-Smirnov 
test applied for the case of Figure lib gives a null-hypothesis prob- 
ability P ~ 3 X 10"". 

One of the clearest separations between the two test samples 
is illustrated in panel (c), where we plot the line shape measure 
kurtosis K.I. in terms of normalized "base" shift. High accretors 
are almost exclusively confined to K.I. < 0.4. The clear excess of 
blueshifted and redshifted bases is evident at low values of K.I. A 
2D Kolmogorov-Smirnov test applied for the two samples gives a 
null-hypothesis P ~ 2 x 10^^^. Panel (d) should be examined in 
conjunction with the previously presented Figure 9c. Figure 12d is 
now restricted to high and low accretors only, as we defined them 
in this section. It is clear now that the two test samples show rather 
different trends on either side of Rpeii ~ 0.50. There is a trend 
extending along the K.I. axis (low accretors) and another one ex- 
panding mostly along the Rfe// axis (high accretors). 

3.1 Composite Spectra of Low and High Accretors 

Starting from Figure 11a we construct median composite spectra 
for two groups of sources. The first group (N=31) is defined by the 
low accreting sources with logLi,o![erg s~^] > 45.5 and normalized 
c(l/4) larger than 0.07 (which is the estimated typical lo uncer- 
tainty of c(l/4)/FW02.5M). The second one (N=26) is drawn with 
the same luminosity constraint from the highly accreting sources 
having the normalized c(l/4) less than -0.07. Focusing on a com- 
mon luminosity regime implies that we are basically emphasizing 
a difference in BH mass. Thus, the composites reflect also the in- 
fluence of the BH masses that are on average different by Idex in 
the two groups (low accretors, mean and median log[MsH(M0)] 
are both 9.3, standard deviation of the mean 0.3); high accretors, 
mean and median log[MflH(M0)] are both 8.2, standard deviation 
of the mean 0.4). We present in Figures 12-13 the composite spec- 
tra for the two cases. We also show our choice for the underlying 



power law continuum and the best solution of the Fell template 
that matches the optical Fell emission, with a focus on the region 
4100-5150A. 

In Figure 12 we analyze the H/J and adjacent emission lines 
in the range ~ 4600-5 150A for the low accretors, after we have 
removed the contaminating Fell in this region and the continuum 
model. We fit simultaneously the lines following the approach 
explained/outlined in earlier in §§ 2.4. We show two possibili- 
ties, first Gauss-l-Gauss for BLR-l-VBLR contributions and second 
Lorentz-nGauss for BLR-l-VBLR. The is very similar in the two 
cases. The BLR component in either case shows a blueshift of ~ 
-500 km s^^ relative to the adopted quasar restframe, based on 
NC. The VBLR on the other hand shows a big difference. In the 
Gauss-HGauss case it has a contribution of ~ 70% to the whole 
Balmer BLR-l-VBLR emission and shows a redshift of ~ -1-2500 
km s~^. In the case Lorentz-l-Gauss the VBLR contributes less to 
the total broad line emission (~ 30%) and is redshifted by ~ -1-4000 
km s~\ The FWHM of the BLR component is close to 5000 km 
s~^ in each case, much less than the FWHM for the combined 
BLR-^VBLR of ~ 7600 km s'^ 

Figure 13 presents the H/3 profile extracted from the com- 
posite of the luminous and highly accreting sources. Because we 
see an inflection along the blue wing of H/3 (indicated in the fig- 
ure) we are attempting a two-component fit (Lorentz-FGauss and 
Gauss-l-Gauss). This time we isolate only H/3 because it is narrow 
enough so that after Fell subtraction its wings go down to the de- 
fined underlying continuum. We also assumed that a narrow com- 
ponent may be present. The solution where the "classical" BLR is 
modeled with a Lorentzian function provides a lower (better by 
a factor 1.6 relative to the Gaussian choice). The blueshift of the 
second component is about -1000 km s~^ and its contributes ~ 
38% to the total broad H/3 emission. In the case of Gauss-l-Gauss 
the second component is blueshifted by -600 km s~^ and has 
a larger contribution^^ 62%. We note that in the fitting procedure 
presented in Figure 13, even though we did not constrain the com- 
ponents to have zero-level baseline, the returned solution is still 
reasonable. 



4 DISCUSSION 

We presented a dual approach in exploring the BLR properties in 
low z quasars: using composite/median spectra and defining a di- 
agnostic set of line profile measures in individual sources. 

The first approach allowed us to reveal that the adopted con- 
cept of two populations A and B is supported by large differences 
in the line profiles of the H/3 emission lines (Figures 3 and 4). Com- 
parison of two representative bins defined in the context of 4DE1 
clearly illustrates that Population B sources typically require an 
additional redshifted component to describe the total H/3 profile, 
while in Population A sources H/3 lines are more symmetric and 
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are best described by a single component. Individual H/3 profiles 
on the other hand show a large diversity displaying shoulders and 
bumps (Figure 5), which are not seen in composite/median spectra. 
This may be telling us that there is a stable emission component 
signature (seen in composites) and a variable one (reflected in the 
large differences between individual and average spectra). 

One of the most important results explored is that low accret- 
ing sources (representative of the so-called Population B) show a 
typical red asymmetric and redshifted Balmer profile at 1/4 frac- 
tional intensity, while the highly accreting sources (found in popu- 
lation A) typically show a blue asymmetric and blueshifted H/3 at 
the profile's base. 

4.1 Very Broad Line Region 

The red asymmetry in Population B can be linked to the proposed 
existence of a distin ct emitting region, the so-called VBLR (e.g., 
iMarziani et ai]|2009h . In that study it is shown that the VBLR con- 
tribution to the total Balmer bro ad emission scales with the source 
luminosity (see also lNetzer & T rakhtenbrot 1Q <S% and becom es in- 
creasingly dominant at higher redshift. Marziani et ai] ( |2009|) also 
suggest that the asymmetry index A.I.(l/4) correlates with the BH 
mass for the intermediate redshift sources (more luminous than the 
SDSS sample). We plot in Figure 14 the asymmetry index A.I.(l/4) 
and the normalized c(l/4) centroid shift versus BH mass (panels 
(a) and (b), respectively), indicating the two populations A and B 
(separated at FWHM(H/3) = 4000 km s"^) with distinct symbols. 
All together (A-l-B) our low z sample does not show a correlation 
between the asymmetry/shift and BH mass. However, there may be 
two distinct trends for the two populations A and B separately. For 
Figure 14a statistical tests Pearson and Spearman indicate for Pop- 
ulation A rp ~ -0.14 (P~ 0.02) and rs ~ -0.17 (P~ 0.006) and for 
Population B rp ~ 0.28 (P~ 3.6 x 10"^), rg ~ 0.29 (P~ 1.9 x 
10"^). For Figure 14b we get for Population A: rp ~ -0.08 (P~ 
0.21) and rs ~ -0.09 (P~ 0.15); for Population B: rp ~ 0.30 (P^ 
8.4 X lO"**), rs ~ 0.35 (P~ 2.3 x 10"^). The statistical tests seem 
to indicate that only Population B sources may show weak corre- 
lations between asymmetry/centroid shift at 1/4 fractional intensity 
and the BH mass. 

If there is indeed a distinct component present in popula- 
tion B sources then the magnitude of its redshift, typically 1000- 
2000 km s~^ is unlikely to be explained by gravitational redshift 
jMarziani et al . 2009, but see also lCorbin 1995). In Fi gure 15 we 
consider the subset of N=249 sources for which we measured a 
positive centroid shift (redshift) at 1/4 fractional intensity (which 
could be interpreted as a sig nature of infalling m aterial). The pre- 
dicted gravitational redshift jPopovic et aljl99^ : 



is much lower than the measured value for virtually all sources. 
Herein Ao = 4862. 7 A is the vacuum wavelength of Balmer H/3, c 
is the vacuum speed of light and G is the gravitational constant and 
the ratio of the BH mass (M) and the BLR radius (R) for a given 
source is estimated as in iSulentic et al.i (i2006) (their formulas 4 and 
7): 

ig. 1.564. 10^3 /™M(H^)y(^^^-r) 

R \ lOOOfcms-i y I \ I 

A pertinent argument in favor of a separate and distinct emit- 
ting region is the existence of clear inflections on the red wings of 
the H/3 profiles in many individual sources (see for example the 



case of B2 1721+34 in ICorbiij|l997h . In Figure 16 we show an 
interesting case, SDSS J230443.47-084108.6, for which we have 
available recorded spectra at two different epochs, about eight years 
apart. The top panel shows the spectr um of the H/3 region observed 
in May 1993 jMarziani et alj|2003ah and the bottom panels show 
the H/3 and Ha spectral regions of the same source, provided by 
SDSS, observed in December of 2001. Not only do we see a sig- 
nificant change in the functional shape of the Balmer profiles, but 
we note the dramatic change on the red side of these lines. In the 
SDSS spectrum we see a clear inflection on the red wing of both 
H/3 and Ha (see also the detailed monito ring study and comp lex 
behavior of NGC 4151 as presented in Sha povalova et alj2009l) . In 
many other cases an inflection between the two components BLR 
and VBLR may or may not clearly show up in the line profiles. The 
two components may be structured differently and apparently they 
do not respond in tandem to the variations in the continuum (see 
the interesting case of PGM 16- 129 in Sulentic et al. 2000^ nor 
do they respond similarly to luminosity variations ( Marzia ni et al.l 
1 2009.) . The search for a clear inflection on the H/3 red wing is fur- 
ther complicated by two additional emission features: i) optical Fell 
has a nontrivial contribution in this region especially due to its mul- 
tiplet m42 components and ii) possible contamination by Hel broad 
permitted emission lines that m ay contribute to the build-up of a 
red-sh elf for both H/3 and Ha (jv eron. Goncalves & Veron-Cettvl 
l2002h . as mentioned earlier. 

However, the presence/absence of inflections is not necessar- 
ily offering us a definitive answer for or against multiple compo- 
nents of the broad lines. The inflections are sensitive to kinematics 
and/or orientation. 

4.2 Disc Wind 

The blue asymmetry/blue-shifted signature in high accretors (Pop- 
ulation A sources) could indicate an additional wind contribution. 
In such highly accreting so urces the presence of r adiation driven 
winds is very likely (e.g., iKing & Poun"d3 12003| ; | Komossa et al.l 
1200 8i) . We note that this potential second contribution on the blue 
wing of H/3 becomes more clearly inflected relative to the clas- 
sical BLR when the [OIII] lines are weak and show a significant 
blueshift, too (see the examples in Figure 5, where the [OIII] shift 
is indicated and see also the more extreme case shown in Figure 
17, where the [OIII] blueshift is ~ 1500 km s~^ and the inflection 
on the blue wing of H/3 is unquestionable). A kinematic linkage be- 
tw een the [OIII] lines a nd the HIL line CIV blueshift was suggested 
bv lZamanov et alj ^2002) . The same kind of linkage may exist be- 
tween the [OIII] lines and the LIL Balmer lines. This idea is further 
reenforced by our Figure 18. We show that there is a correlation 
between the shift of CIV (rest frame AI549A) at 1/2 fractional in- 
tensity and the normalized shif t of H/3 at 1/4 fractional intensity 
(see also similar comparisons in lMarziani et al.|[l996) . The data in 
the plot represent the common set of sources of the current sam- 
ple of SD SS quasars and the sample explored spectroscopically in 
Sulentic et al. (2007). Although the size of the overlapping sample 
is rather small (n=25) there is an important observation: there is no 
source showing a CIV redshift and an H/3 blueshif t. We also note 
that a c omparison between the Atlas presented in IMarziani et al. 
("2003 a) for the H/3 and the sample presented in ISulentic et al.l 
(2007) for CIV reveals an overlap of n=87 sources. 69 of those 
sources (~ 80%) show a shift of CIV that is larger than the shift in 
the Balmer line (at both 1/4 and 1/2 fractional intensity). This kind 
of empirical evidence suggests that the blue asymmetry/shift of H/3 
might be produced by the same mechanism that is responsible for 
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CIV shift, and the simp l est interpretation may be a wind/out flow 
(e.g. jMurrav et alJll995l : lKollatschnvll2003l:lYoung et al.l2007h . 

The shift of quasar broad emission lines (relative to their rest- 
fame) has_been_jntensively studie d (e.g., iGaskell 1982; Wilke: 



1984 ICarswell et al 



^_^[T991; Tvtle r&Fanl TT992: Corb inI 1 19951: 
Mcintosh et al.l I199S1 : IVanden Berk et al.l I2OOI. : ,Richards et al.l 
2002h . It has been reported that the amplit ude of the shift is cor- 



1991; Tvtler & FanI 



1992; 



related with the ionizati on potential (e.g., lMcIntoshetai]|l999l ; 
IVanden Berk et all 1200 ih . The aforementioned comparison be- 
tween the shift of the H/3 and CIV (based on the Atlas of 
Marziani et alj20033 and the sample of quasars from lSulentic et all 



2007h 



supports it. 



4.2.1 Other Arguments for the Wind Scenario 

iKallman & Kroii^ ( Il986l) argued that when HILs are systemati- 
cally blueshifted relative to the quasar's restframe by thousands 
of km s~^ it becomes more likely that absorption lines with Zabs 
> Zem will exist. Let us examine the concrete case of quasar PG 
1700-1-518, which shows both low and high ionizat ion absorptions 
jLaor & Brand3 12002| ; [Pettini & Boksenberg 1985) a nd a clearly 
blue-a syrmnetric H/3 profile (^arziani et al. 2003a). lYoung et al.l 
feOOTl) offer arguments in favor of a strong rotating wind in this 
source. This scattering wind rises nearly vertically from the disk 
over a range of radii approximately c oincident with the B aimer 
emitting region. So, the assertion of Kallm an & Krolilj l ll986l) 
could be extended to the blueshifted components of LILs like 
H/?, i.e., when sources show large blue-asymmetries in the Balmer 
emission lines, the BAL signatures in Mgll A2800A or CIV 
AI549A could be present as well. In the case of PG 1700-1-518 it 
is also possibl e that we see th e BLR disk more face-on than equa- 
torially ( Richards et al.ll2002h and the blue-asymmetry is formed 
by scattering in a wind out flow. Another similar case studied in de - 
tail is IRAS 07598-H6508 jLiparilll994 IVeron-Cettv et al.ll2006l) . 
This source shows both Mgll and CIV BAL signatures and an in- 
flected blue wing of the H/9 emission line. Both quasars mentioned 
here have been interpreted as examples of nuclear young/starburst 
(IRAS 07598-1-6508) or poststarburst (PG 1700-H518) by iLiparil 
( I1994I) . 

Exploiting the analogy between microquasars, stellar-scale 
phenomenal and and quasars, galactic-scale phenomena, one can 
also interpret the recent studv lNeilsen&Led 120091) as supportive 
for the idea that the presence of disk winds inhibits the formation 
of jetted radio structures. This may have implications for the fact 
that FRII quasars and sources where disk wind signatures are more 
evident occupy mutually exclusive regions in the optical plane of 
the 4DE1 parameter space. 

4.3 Other Scenarios for the Shifted Broad Lines 



[Richards et ^ j2002l) offer a different view on the matter of the 
CIV blueshift. They propose that the blueshift is an apparent effect 
produced by obscuration of the redshifted photons. In the proposed 
scenario the largest blueshifts (which may not be real, but rather 
due to the lack of flux in the attenuated red wing) are explained by 
an equator ial viewing angle of t he line emitting region(s). An ar- 
gument of [Richards et alj j2002l) that favors an equatorial view of 
the BLR when we observe increasingly larger blueshifts of the HIL 
like CIV is that the composite spectrum of the most blueshifted 
CIV spectra is most similar to that of BAL quasars. This similar- 
ity would indicate that the red wing of CIV is attenuated by ob- 
scuring material. Another argument raised in favor of obscuration 



coupled with orientation (external or internal) is the finding that the 
strength (equivalent width) of the CIV line declines with increasing 
blueshift. This is difficult to reconcile with the Fell emission that 
is likely produced in a high density/column density environmen t 
(flattened geometry; e.g.. lCollin-Souffrin. Hameurv & Jolvll988ah . 
We do observe that the average strength of the H/3 broad emission 
line declines by a factor two from Al to A4 (Figures Ib-c) and 
the H/9 blue asymmetric and blueshifted base is typical for bins 
A3/A4 (lower equivalent width W(H/3)), however, this can be ex- 
plained by a quenchin g of the Balm er emission in high electron 
density environments jGaskellllioOol) , which favor the Fell emis- 
sion (as we see in Figures Ib-c, W(Fell) increases by the same 
factor two from bin Al through bin A4). Largest blueshifts cor- 
respond to naiTowest and strongest Fell, and we suggested that 
Fell decreases from f ace-on because of a sort of "limb darkening" 
jMarziani et al.|[l993) , which supports a wind scenario plus disk 
obscuration, but a different orienta tion dependence from the one 
suggested by [Richards et alj j2002l) : the largest blueshifts are ob- 
served when the source is viewed pole-on. 

[Richards et"!!] ( [2002[) advocate the ubiquity of the CIV 
bl ueshift. This la st result is not confirmed by the most recent study 
of [Sulentic et alj (2007), where only Population A quasars show 
such a spectral signature. Population B sources (where the major- 
ity of RL quasars belongs) show unshifted and sym metric CIV pro- 
files. Most likely the bins of least CIV blueshift in [Richards et al.[ 
(12002), which also have the largest radio-detected fraction and the 
largest CIV equivalent width (see their Table 2), must have a sig- 
nificant overlap with the Population B sources. It is difficult to see 
how attenuation of the CI V red wing, ac c ording to the orientation 
dependence suggested by [Richards et alj ( [2002I) . can leave a sym- 
metric profile in so many quasar s: population B sources account for 
- 60% of our CIV sample tSulentic et al.[[2007[) . 

An alternative to the the wind interpreta ti on of blueshifts was 
recently proposed by [Gaskell & Goosmand j2008l) . They advo- 
cate that the blueshift of HIL like CIV is produced by Rayleigh 
and/or electron scattering off of inflowing material of the BLR. 
Their model predicts that the relative blueshifts of different lines 
in the same AGN should be proportional to the inverse square 
root of the radius where they are produced. The results we indi- 
cated in the previous paragraph, that the H/3 blueshift is almost 
al ways less than that of CIV are consistent with the prediction 
of [Gaskell & GoosmannI ^20081) in the sense that CIV is proba- 
bly emitted at a smaller radius within the BLR relative to the 
Balmer emitting region. This is qualitatively supported also by the 
fact that HIL CIV line shows broader FW HM profiles than H fl 
in Popula tion A sources (but not in B ; see [Sulentic et alj[2007h . 
Moreover ICaskell & GoosmanrJ ( I2OO8I) predict the amplitude of 
the blueshift is proportional to the infall velocity. Spectroscopic 
signatures of the inflow ha ve been indicated (see section 2.1 of 
[Gaskell & Goosmand 20081. see also .Bentz et al. 2008), although 
very recent reports jPermev et al.[[2009l : Ishapovalova et al. | [2009l) 
indicate evidence for a very diverse behavior of the BLR, which can 
show signatures of outflowing, inflowing and virialized gas motion 
in different sources. 

If the scattering effect depends on the wavelength only, then 
low-ionization lines in UV should show a blueshift c omparable to 
that of the UV high-ionization lines. Published data [Tvtler & FanI 
(l992) as well as a tentative analysis on several HST spectra (Ne- 
grete et al. in preparation) do not support this suggestion. The 
01 1304 A for example, does not seem to show the large blueshift 
of CIV and seems to be more consistent with the quasar rest-frame 
as derived from the narrow low-ionization emission lines. 
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Recently IHu et iD l l2008allbh propose that inflow may be a 
ubiquitous phenomenon as reveale d by their measur ed (red-) shift 
of the Fell optical lines. However, I Hu et alj ( l2008al lbl) report that 
the infall is increasingly pronounced as the accretion rate (or Ed- 
di ngton ratio) decreases. This w orks against the model proposed 
bv lGaskell & GoosmannI | |2008|) because the largest blueshifts are 
typically observed in highly accreting sources (members of the 
Population A). The lHu et al.l ( l2008al lbl) recent papers propose that 
there is an Intermediate Line Region (ILR) component in H/3 that 
is systematically redshifted and narrower than a so-called VBLR 
rest-frame component. This is a rather different view compared 
to the proposed two-c omponent model of ISulentic et alj ([2002) or 
iMarziani et al. I ( |200^ . where the VBLR shows a significant and 
systematic redshift and the classical BLR signature is basically 
at the quasar's rest frame. A two-compon ent model employin g 
the VBLR/ILR terminology is proposed by J Popovic et ani2004l) . 
where the emission of the line wings comes from a disk and the 
line core from a distinct emitting region that has spherical geome- 
try. The two emitting regions may be linked by a wind. One must 
be cautious about different definitions and labels of the broad line 
emitting regions used by different studies; the same name may be 
assigned to different structural components. 

We are not y et convin c ed of t he ubiquity of the Fell red- 
shift proposed by IHu et alj feOOSdlb h. Although we don't dis- 
miss that in sources like SDSS J094603.94+013923.6 or SDSS 
J101912.57-I-635802.7 (both shown in our Figure 5) the Fell could 
be redshifted, we also encounter sources like SDSS J001224.02- 
102226.£] or SDSS 1094756.00-^535000.3 that can be difficult to 
fit in the scenario they propose. And we also find an excess of 
blueshifted peaks in both populations A and B, which does not sup- 
port the idea of a redshifted ILR in virtually all quasars. This is fur- 
ther confirmed by the composite spectrum shown in our Figure 4 
for bin B2, which should show a significantly redshifted Fell signa- 
ture, which we do not see. Moreover, we selected the (apparently) 
brightest sources in their sample (psf i or psf g <17.5, just as we 
const rained our SDSS low z sample; IZamfir. Sulentic & Marziana 
l2008i) . We found about 50 such sources with a claimed Fell red- 
shift larger than 1000 km s~^. The authors state that their sources 
are bright enough that the host galaxy plays a negligible role. 
We did find sources where prominent absorption lines indicate a 
host galaxy contamination, e.g., SDSS Jl 10538.99-1-020257.3. This 
raises the question whether the mismatch between the Fell rest- 
frame template and the spectral lines is severely affected by the 
dips created by an ignored host galaxy signature. The problem 
has profound implications toward understanding the link between 
the 4DE1 space (through RfeJ/) and i ts proposed main d river 
Lboi/Lfidd (see the very recent study of iFerland et al.ll2009l) and 
thus requires further scrutiny. 



4.4 The Relative Fraction of Population A/B 

IShen et alj bOOSl) compile BH masses for 60,000-1- Sloan 
quasars within redshift z =4.5, starting from the catalog of 
ISchneider. D.P et aljJioOTi) . which imposes a faint limit of -22 for 
the absolute i-band mag nitude. In or der to compare to our sample, 
we selected from iShen et al.l ( l2008h (using the VizieR catalogue 



^ A similar case SPSS J 153636.22-I-044127.0 was found in SDSS DRV by 
iBoroson & LaueJ j2009 b and interpreted as a possible signature of a sub- 
parsec binary BH system, although a number of other articles propose al- 
ternative explanations. 



database, lOchsenbein. Bauer & MarcoutI |2000|) only the quasars 
with z<0.7. There are 12,003 listed sources, of which 8,979 have 
reported measurements of the H/3 line. We were unable to under- 
stand why ~ 25% of the sample has no reported measures of H/?. 
They are neither fainter, nor noisier/more problematic spectra than 
the rest. The other major criticism is the relative proportion of 
sources that we would call A/B in the light of our adopted 4000 
km s~^ boundary. They report 71% sources with FWHM(II/3) > 
4000 km s^^ (of those 8,979 measured). We tested whether and 
how this proportion is tied to the S/N reported by their study. We 
found that as the S/N decreases there is an increased fraction of 
Population A sources (i.e., 19% for S/N>20, 26% for S/N>15- 
20, 30% for S/N>5-15). All previous m ajor studies of quasars in 
the context of Eigenvector 1 space (e.g., B oroson & GreenI 1992 ; 



Sulentic. Marziani & Dultzin- Hacvanl2000al ; ISulentic et al.l2000bl; 
Marziani et alj2003ah report a slight numerical dominance of Pop- 



ulation A sources (at low z), i.e., FWHM(H/3) < 4000 km s"\ Our 
main concern is their adopted H/3 fitting routine, i.e., a single Gaus- 
sian for all profiles. We now know that Population A spectra are 
best described by Lorentz-profiles, while for Population B sources 
is a double-Gaussian model works better, due to their pronounced 
re d asymmetry ^S ulentic et al. 2002; Netzer & Trakhtenbrot 200^; 
Mar ziani et alj2009l) . A single Gaussian fit would overestimate the 
FWHM measure in both Population A and B sources; for instance 
a single Gauss fit to the A2 or Bl bin median composites spec- 
tra (discussed earlier in §§ 2.4) would overestimate the FWHM of 
the broad Hbeta by 20-25 %. Another serious problem with their 
analysis is the inclusion of a large number of spectra of sources 
that are either in a high-continuum phase are the emission lines are 
intrinsically very weak (e.g., SDSS J083353. 88-1-422401. 8, SDSS 
J083 148.87-1-042939.0), where the broad lines are almost miss- 
ing. For such cases they report hugely overestimated values of the 
FWHM(H/3). 



5 CONCLUSIONS 

We have shown that the Population A/B concept proposed in the 
context of the 4DE1 Parameter Space holds robustly when other 
line diagnostic measures like centroid shifts, asymmetry and kur- 
tosis are considered. Taken at face value the spectral measures tell 
us that populations A/B show a great deal of difference in terms 
of BH mass and Eddington ratios (Figure 10). This may be indi- 
cating that the numerous observed dissimilarities between A and B 
are indeed driven by large differences in the structure and kinemat- 
ics of the central engine. Nonetheless, large differences in terms 
of BH mass and Eddington ratio may be reduced. On one hand it 
was pointed out that NLSyl sources may not harbor such small BH 
masses as inferred from their narrow emission lines. They might 
not be super- Eddington ac cretors if proper correct ion for radia - 
tion pressur e jMarconi etal . 2008, 20 0^, but see alsolNetzei<2009l) . 
orientation ("Marziani et al. 2003h; iDecarli et al.ll2008l) and geom- 
etry (e.g., Collin et al. 2006) are made. On the other hand, if the 
H/3 profiles that show a large red asymmetry, which is probably 
caused by a distinct VBLR redshifted contribution, are properly 
corrected for such "unvirialized" contributions, then the estimated 
BH mass gets smaller and thus the measured Eddington ratio gets 
higher. The large variance of the 4DE1 space would be signifi- 
cantly reduced from both sides of the accretion rate sequence. As 
a consequence, the FWHM measure, which is an important op- 
tical dimension within 4DE1, may be replaced with A.I. (1/4) or 
the normalized c(l/4)/FW025M centroid shift. They correlate at 
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least equally well with Rfe/j (Figure 9). Nonetheless, incorporat- 
ing these two new measures in the 4DE1 is not an easy task, as they 
also manifest ( at least in Population B) a sensitivity to luminosity 
jNetzer & Trakh tenbrot 2007; Marziani et al. 2009 ), which was in - 
dicated as part of the Eigenvector 2 space feoroson & Greejl992h . 
formally orthogonal to Eigenvector 1 . Moreover, the Population B 
H/3 asymmetry/centroid shift at the base (1/4 fractional intensity) 
appe ars weakl y corr elated with the BH mass (Figure 14; see also 
iMarziani et aU2009l) . 

We also presented several arguments that two populations 
of quasars may be redefined based on an Rfe/i nominal bound- 
ary of 0.50. This could be seen as an alternative empirical defi- 
nition of some equivalent Populations A and B. The lines of ev- 
idence are as follows: 1) RFe//=0.50 seems a critical value that 
best iso lates RL quasars from the rest of quasars (recall Figure 
4 from 'Zamfi r. Sulentic & Marzianil 1200^ . 2) around this mea- 
sured value of Kpeii the centroid shift and the asymmetry of H/3 
change sign (Figures 9a-b), 3) a plot of the line shape kurtosis in- 
dex K.I. versus Rfe// shows two different trends on either side 
of RfeJ/ ~ 0.50 (Figure 9c), likely separating low from high ac- 
creting sources (Figure 12d) and 4) a tight correlation between 
Kpeii and the Eddington ratio indicates that Rfe// ~ 0.50 roughly 
corresponds to Lboi/LEdd ~ 0-18 (Figure 9d), consistent with 
several other studies that suggest a critical change in the struc- 
ture/geometry/kinematics in the quasar "engine" around that value. 
The red and blue H/3 base asymmetries observed at the opposite 
ends of the Eddington ratio sequence may have very different ori- 
gins. 

We also find that our measures for low accretors (typically 
members of Population B) seem to match the predictions of accre- 
tion disk models, i.e., blueshifted peaks c(9/10) and redshifted/red 
asymmetric bases c(l/4) and A.I.(l/4) (Figure 12b). This allows 
us to suggest a typical flattened BLR geometry for all Population 
B quasars, not only for the RL quasars, whose majority is seen in 
Population B subspace. We emphasized that our conclusions are not 
driven by radio-loudness. Focusing our investigation on the radio- 
quiet subsample the results are similar to what we get using the 
whole sample of N=469 quasars. 

Although the sample of FRII sources is relatively small 
(N=44), we find a very interesting correlation between the 
core:lobe flux density ratio, which may be seen as an orientation 
indicator (based on the 1.4GHz FIRST survey measures) and the 
redshift of the H/3 base (Figure 19). The more core-dominated FRII 
sources show a more pronounced redshift at 1/4 fractional inten- 
sity (see also Figures 9a-b and related commentaries). We may just 
speculate that the excess of redshifted photons could be a signa- 
ture of an inflow feeding the central mass and allowing for a jetted 
bursB A larger sample should be investigated. 
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Figure 1. (a) The optical plane of the 4D Eigenvector 1 Parameter Space. The quadrants indicate the bins defined bv lSulentic et alji2002l) and adopted here for 
composite spectra. Typical 2cr eiTor bar s are shown in the upper right comer. The horizontal dash-dot line indicates the nominal boundary between Populations 
A and B. The relative numerical contribution of from each population is 55% in A and 45% in B. (b) The number of sources considered in each bin along 
with a few (median) relevant measures in each bin are indicated. The square in the upper light is the legend of the numbers shown in each bin. Inside each bin 
(lower-right square) we report the number of FRII and core-dominated RL sources along with the estimated radio-loud fraction, (c) Similar to (b), this time 
the numbers represent mean and standard deviations (in the parenthesis in each case). 




Figure 2. Distribution of the S/N measure (per pixel) in the 5600-5800 A region for our sample of N=469 sources. 
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Figure 3. (a) Median composite spectrum of A2 bin sources and tlie best model of Fell emission around H/3. (b) A simultaneous fit (our best solution) of H/3 
and adjacent lines in the selected wavelength interval of the composite. 




Figure 4. (a) Median composite spectrum of B 1 bin sources and tlie best model of Fell emission around H/3. (b) A simultaneous fit (our best solution) of H/3 
and adjacent lines in the selected wavelength interval of the composite. 
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Figure 5. A selection of H/3 profiles meant to illustrate the diversity of asymmetries, shapes and shifts in individual sources. A magenta dash-dot vertical 
line in each case indicates the rest-frame defined by H/3 NC at A4862.7A (Note that we refer our measures to vacuum rest-frame wavelength values). A blue 
dotted vertical line indicates the peak of the broad H/3 defined at 9/10 fractional intensity. In a few panels a third vertical line indicates the rest-frame [OIII] 
A5008.2A and the actual blueshifted position of that line. The displayed spectra still include the Fell contribution, but the H/3 measures employed in the paper 
are performed on FeII-"cleaned" profiles. Panels are grouped by bins and displayed in order of increasing asymmetry index A. I. (1/4) within each bin. In each 
case we also provide a set of numbers: bin, bolometric luminosity, S/N measures in the continuum (see Figure 2), asymmetry and kurtosis indices, centroid 
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Figure 5. Cont'd. 
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Figure 6. Plots of FWHM(H/3) versus (a) line asymmetry, (b) shape measure kurtosis, (c) "base" normalized shift c(l/4)/FW0.25M and (d) "peak" normalized 
shift c(9/10)/FW0.90M. The vertical lines indicate the 2cr uncertainty on either side of zero in all but panel (b). In case we also show the median typical 2ct 
error bars. The horizontal line at 4000 km s~ is the Population A/B boundary proposed in the context of the 4DE1 space. The FRII sources are shown with 
distinct blue symbols in each case. 
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Figure 7. Histogram distributions of the asymmetry index A. I. (1/4) - upper left, base centroid shift c(l/4) normalized to FW0.25M - upper right, peak centroid 
shift c(9/10) normalized to FW0.90M - lower left and line shape kurtosis index K.I. - lower right. We show separately the distributions for Population A and B 
sources for each measure. The width of the bins is equal to a estimated typical eiTor. The vertical lines indicate a ±2(T interval on either side of for the first 
three parameters. 
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Figure 8. Histogram distributions for tlie centroid shifts at the base (1/4 fractional intensity) and peak (9/10 fractional intensity) without normaUzation to the 
width of the profile. Population A and B are illustrated separately. The width of the bins is equal to a estimated typical eiTor The vertical lines indicate a ±2(t 
interval on either side of in each case. 
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Figure 9. The connection between the hne diagnostic measures and the Eigenvector 1 optical parameters and its chief driver, Eddington ratio. Darker symbols 
indicate (in each panel) the FRII sources of our sample, assumed to be the parent population of RL quasars along with typical 2rj error bars in the comer of 
the panel, (a) The asymmetry index is plotted as a function of Rj?^//. The bisector fit is also indicated, (b) The normalized base centroid shift is plotted as a 
function of Rfe/i ^nd the bisector fit is also shown, (c) The line shape measure kurtosis index K.I. is plotted as against R^ejj. The vertical line suggests a 
possible separation of two different trends around 0.50 abscissa value, (d) Eddington ratio is shown as a function of RfcJI along with a bisector fit. Panels 
e-f show the relation between base measures A.I.(l/4) and normalized centroid shift c(l/4)/FW0.25M and the Eddington ratio, the main driver of the 4DE1 
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Figure 10. Population A and B quasars show similar luminosity distributions, but strikingly different distributions in terms of BH mass and Eddington ratios. 
Population A typically showing smaller BH masses and higher accretion rates and population B harboring more massive BH and accreting much slower 
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Figure 11. Two subsamples of N=177 low (X^boll^Edd < 0-1) andN=I23 highly accreting Q^bol/^Edd > 0-3) sources are considered, the squares representing 
the former and the circles showing the latter. The filled symbols indicate the FRII sources. Pairs of vertical and horizontal lines on either side of zero show the 
typical ±2(7 median uncertainty, (a) The bolometric luminosity versus the base (1/4 fractional intensity) normalized shift. The horizontal line at 45.5 ordinate 
indicates the values above which both samples are well represented numerically (see text), (b) The normalized peak shift versus normalized base centroid 
shift, (c) Kurtosis index versus normalized base centroid shift. Typical median 2a errors are shown in the upper comer of the panel, (d) Kurtosis index versus 
Rpell- Typical median 2cr errors are shown in the upper comer of the panel. 
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Figure 12. Upper-panel: Composite median spectrum for N=31 sources with logL(,o;[erg s~^] > 45.5 and c(l/4)/FW0.25M > 0.07. Tlie best model of Fell 
is also shown along with the underlying continuum. Lower-left and lower-right panels show the best fitting solution for the emission lines in the interval 4600- 
5 150 A. The left panel shows a Gauss-Gauss combination for the H/3 profile, the right panel shows a Lorentz-Gauss model of H/3. All other lines/components 
are modeled with Gaussian functions. In both panels we also show the residuals. 
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Figure 13. Upper-panel: Composite median spectrum for N=26 sources witli logL5o;[erg s^^] > 45.5 and c(l/4)/FW0.25M < - 0.07. Tiie best model of 
Fell is also shown along with the underlying continuum. Lower-left and lower-right panels show the best fitting solution for the emission lines in the interval 
4600-5150 A. The left panel shows a Gauss-Gauss combination for the H/3 profile, the right panel shows a Lorentz-Gauss model of H/3. In both panels we 
also show the residuals. 
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Figure 14. Population A and B sources (open circles and solid squares, respectively) are shown in plots of BH mass versus asymmetry index A. I. (1/4) (panel 
a) and normalized centroid shift c(I/4)/FW0.25M (panel b). In each panel we also indicate typical 2cr en'or bars. 
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Figure 15. For N=249 sources that show a positive centroid shift c(l/4) we plot that measured shift (solid squares) and the predicted gravitational redshift 
(open circles) versus the FW0.25M. We also indicate the typical 2ct error bars for the measured values of c(l/4). 
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Figure 16. SDSS J230443.47-084108.6 observed at two different epochs. We note the change in the functional shape of the Balmer H/3 hne wings. For the 
more recent records (lower panel) we show both the H/3 and the Ha regions. The vertical lines indicate the rest-frame Balmer wavelength positions. The 
"redshelf" is clear in the later observation, especially in the Ha region. An inflection along the red wing of the Balmer lines may be a transient feature due to 
the variability of the emitting source. 




Figure 17. An extreme case that shows a clear inflection on the blue wing of H/3 accompanied by a -1500 km s ^ blueshift of [Olll] A5008A. The vertical 
lines indicate the restframe positions for H/3 4862. 7A and [OIII] 5008. 2A. 
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Figure 18. A compa rison between the CI V A1549A and Balmer H/3 shift. Data shown represents the overlap between the current sample and the sample of 
quasars examined in lSulentic et alj i2007h . The dotted lines indicate the zero shift for CIV (vertical axis) and H/3 (horizontal axis). Typical la uncertainties 
are indicated. 
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Figure 19. The FRII sample is shown here in terms of core:lobes flux density ratio at 1.4GHz, based on FIRST and the normalized H/3 emission line base 
shift. The vertical line indicates an unshifted measure. 



